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We investigate the transition from unitary to dissipative dynamics in the relativistic 0{N) vector 
model with the A(</9^)^ interaction using the nonperturbative functional renormalization group in 
the real-time formalism. In thermal equilibrium, the theory is characterized by two scales, the 
interaction range for coherent scattering of particles and the mean free path determined by the rate 
of incoherent collisions with excitations in the thermal medium. Their competition determines the 
renormalization group flow and the effective dynamics of the model. Here we quantify the dynamic 
properties of the model in terms of the scale-dependent dynamic critical exponent z in the limit of 
large temperatures and in 2 < d < 4 spatial dimensions. We contrast our results to the behavior 
expected at vanishing temperature and address the question of the appropriate dynamic universality 
class for the given microscopic theory. 
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I. INTRODUCTION 

Efforts to understand the properties of extreme states 
of matter in QCD mm rely crucially on the identifica¬ 
tion and characterization of universal fluctuations. Un¬ 
derstanding their properties is essential to settle the pres¬ 
ence and nature of the critical end point [3] , marking the 
end of the first order transition between the hadronic and 
quark-gluon plasma (QGP) phase. To establish the hy¬ 
drodynamics of the QGP near the QCD critical point, 
it is in particular their dynamic properties that are of 
considerable interest [1]. Near the critical point the char¬ 
acteristic fluctuations display the phenomenon of critical 
slowing down which has important ramifications: It con¬ 
trols the maximal achievable correlation length in the 
process of the rapid cooling of the quark-gluon plasma 
through the critical region [5] and thereby sets the typ¬ 
ical strength of event-by-event fluctuations in heavy ion 
collisions probed at RHIC or the LHC [5]. On the other 
hand, it constrains important nonequilibrium effects as, 
e.g., the production of topological defects mi- It is 
a remarkable fact that these phenomena are determined 
only by a few quantities characteristic of the universality 
class of the critical point. While the identification of the 
static universality class is dictated by principles of sym¬ 
metry and the dimensionality of the problem mm , to 
identify the proper dynamic universality class requires an 
understanding of relevant modes at the phase transition 
(i.e., the conservation laws respected by the dynamics). 
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Unfortunately, there are no general principles that al¬ 
low to distinguish whether a particular mode is relevant 
at the critical point, or if its effect is irrelevant for the 
dynamic critical behavior. To answer this question it 
is necessary to reconcile well-established effective models 
for critical dynamics used in the vicinity of classical phase 
transitions mm with our understanding of the under¬ 
lying coherently propagating excitations. In the past, 
there have been several attempts in this direction using 
perturbative methods as, e.g., the dynamic renormaliza¬ 
tion group [13 EH, as well as nonperturbative methods, 
e.g., AdS/CFT duality [IKHTT] . classical-statistical lat¬ 
tice field theory simulations [18] , or n-particle irreducible 
(nPI) effective actions [in]- We propose a new attempt 
to make a solid connection between dynamic critical phe¬ 
nomena and microscopic quantum physics using the non¬ 
perturbative functional renormalization group (RG) [20j . 
This formalism has the potential to shed new light on 
the smooth quantum-to-classical crossover and enables 
us to address the question of the appropriate dynamic 
universality class for a given microscopic theory. In par¬ 
ticular, we provide a continuous connection between the 
microscopic relativistic and nonrelativistic dissipative dy¬ 
namics at nonvanishing temperature - going beyond con¬ 
ventional approaches to dynamic critical behavior based 
on effective models. Our work constitutes a first impor¬ 
tant step towards the goal to fully unravel the dynamic 
properties in both the quantum and classical domain in 
the framework of the real-time functional renormaliza¬ 
tion group. 

In this paper, we consider the relativistic 0{N) vector 
model with the interaction and investigate the 

transition from the microscopic dynamics to low-energy 
effective dynamics in the vicinity of a continuous phase 
transition. Let us briefly explain qualitative properties 
of this model. It can be seen to describe an interacting 
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FIG. 1. Illustration of the expected phase diagram for the 
0{N) vector model in d > 2 spatial dimensions. It fea¬ 
tures a line of second-order phase transitions T^r = Ter {5 k) 
(blue, continuous line) that separates an ordered from a dis¬ 
ordered phase, and terminates at a quantum critical point 
(T = (5 k = 0). As this critical line is approached, the corre¬ 
lation length ^ = ^(T, 5k) diverges. Shaded areas highlight 
those regions in the phase diagram where critical fluctuations 
become important. The narrow region around 5k is controlled 
primarily by the QCP and becomes smaller as the spatial di¬ 
mension d is increased. It vanishes completely above d = 3, 
in which case the effect of fluctuations can be neglected. The 
critical region around the thermal transition TcriSn) similarly 
decreases in size until fluctuations cease to be important above 
its upper critical dimension d = 4. 


gas of relativistic bosonic particles in thermal equilibrium 
where the interaction is short-ranged and determined by 
the s-wave scattering length Ultraviolet renormal¬ 
ization tells us that a finite number of parameters, which 
determines the low-energy scattering of particles in the 
vacuum, totally characterizes the microscopic dynamics. 
On the other hand, the screening of soft momenta at non¬ 
vanishing temperature is achieved by infrared renormal¬ 
ization. Thus, there are two competing processes that 
play a role in the system: The coherent scattering of 
particles leads to a renormalization of the coupling and 
mass parameters of the model, whereas incoherent col¬ 
lisions with thermal excitations in the medium endow 
the particles with a finite lifetime. The importance of 
these two processes for the critical behavior of the the¬ 
ory strongly depend on the spatial dimension as well as 
the temperature of the system. 

At zero temperature (T = 0) and in spatial dimension 
d > 2, the relativistic 0{N) model features a second- 
order phase transition which separates a disordered (sym¬ 
metric) phase from an ordered phase with a nonvanishing 
field expectation value. In the vicinity of such a quan¬ 
tum critical point (QCP) [52] the system is character¬ 
ized by a diverging correlation length ^ ~ , where 

(5k = K — Kc is a microscopic parameter that controls the 
deviation from the critical point and v is the correlation- 
length exponent. The diverging correlation length can 
be related to the vanishing of a relevant spectral gap 


A ~ where z defines the associated dynamic critical 
exponent!^ There is no characteristic scale at the QCP - 
the system fluctuates at all scales and thus some physical 
quantities, such as scaling exponents, amplitude ratios, 
etc., can be universally classified without knowing the 
details of the microscopic physics. 

Temperature is a relevant perturbation at the QCP. 
It yields a different universality class for the thermal 
transition or may even render the transition unstable. 
While for d > 2 the zero-temperature transition lies in 
the (d-b l)-dimensional universality class of the classical 
0{N) model, at nonvanishing temperature the system 
experiences an effective dimensional reduction and the 
appropriate characterization of universal fluctuations is 
given by the d-dimensional theory. For spatial dimension 
d = 2 the QCP is in the three-dimensional 0{N) univer¬ 
sality class, however for T ^ 0 one expects that the field 
expectation value always vanishes, in accordance with 
the Mermin-Wagner theorem [55] Indeed, for vector 
models with N > 3 this is the case and there are no 
indications for an ordered low-temperature phase. Nev¬ 
ertheless, for A = 2 the system experiences a Kosterlitz- 
Thouless transition at some critical temperature, which 
separates a high-temperature disordered phase from a 
low-temperature phase with algebraic order. This obser¬ 
vation does not contradict the Mermin-Wagner theorem, 
which cannot exclude a low-temperature phase with al¬ 
gebraically decaying correlations (see, e.g.. Ref. [24] and 
references therein). In the presence of a nonvanishing 
temperature we may identify different regimes for the 
dynamic correlations and response of the system: In the 
classical regime A ^ T, all modes are highly occupied 
up to the scale set by the temperature. On the other 
hand, when the temperature is sufficiently small, the 
thermal occupation of modes is exponentially suppressed, 
i.e., ^ 1, for io{p) > A. This defines the quan¬ 

tum regime, where occupation numbers are small and the 
commutator of field operators cannot be neglected. Close 
to the QCP, in the quantum regime, the dynamic critical 
exponent z = 1 in the relativistic 0{N) model which is 
a consequence of Lorentz symmetry. As the temperature 
is increased and the system crosses over to the classical 
regime, similar to the change in static universality class, 
the dynamic critical exponent changes (z > 1). It is this 
interplay of temperature and dynamics that we address 
in the following. These arguments are summarized in 


^ We only indicate the leading divergence in ^ and A and thus 
the provided scaling holds up to an appropriate dimensionful 
constant (and additional subleading scaling corrections). Note, 
that both the correlation length and the spectral gap depend on 
the temperature T, the parameter 5k, and possibly other external 
fields, in general. 

^ There is an important exception to the theorem and that is the 
theory with a discrete order parameter (TV = 1). At T yf 0 and in 
spatial dimensions d > 1, this model features a continuous phase 
transition in the d-dimensional Ising universality class, while the 
QCP lies in the (d -|- l)-dimensional Ising universality class. 
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Fig. which shows an illustration of the expected phase II. REAL-TIME n-POINT FUNCTIONS AND 
diagram at zero and nonvanishing temperatures. GENERATING FUNGTIONALS 


We provide a picture of the processes that lead to the 
demise of coherent quasiparticle excitations and demon¬ 
strate how an effective classical description of many-body 
systems emerges at nonvanishing temperature. For that 
purpose we employ the nonperturbative functional RG 
m that allows us to study the impact of temperature 
and dimensionality independent of the assumption of 
small coupling or small deviations from the upper critical 
dimension. Similar work along these lines has been pur¬ 
sued, where the properties of classical phase transitions 
and its static critical exponents have been examined in 
the real-time formalism, e.g., at the example of the real 
scalar theory [2S] and gauge theories m Dynamic 
scaling properties have been addressed in the framework 
of the functional RG in the context of effective models for 
critical dynamics |12] . This includes in particular. Model 
A [2H] and Model C [H] in the Halperin-Hohenberg clas¬ 
sification of dynamic universality classes m- Further 
applications include stochastically driven systems out of 
equilibrium [30H32], e.g., reaction-diffusion and percola¬ 
tion problems [33], the Kardar-Parisi-Zhang (KPZ) equa¬ 
tion [33], hydrodynamic turbulence [331|33], and nonequi¬ 
librium steady states in closed systems [371 EH]. Close to 
a quantum phase transition (QPT) the functional renor¬ 
malization group has been applied to a number of prob¬ 
lems in the imaginary-time formalism. These include 
nonrelativistic systems where a crossover was observed to 
a quasi-relativistic action with second (imaginary) time 
derivatives [33] and the two-dimensional 0{N) model [30] 
for which analytic continuation was used to derive spec¬ 
tral functions in the quantum critical regime. The rel¬ 
ative importance of quantum and thermal fluctuations 
in the vicinity of a QPT was addressed for the two- 
dimensional quantum Ising model |41j . as well systems 
with fermions |42j . 


The outline of this paper is as follows: In Sec. [TT]we re¬ 
view real-time formalisms to determine n-point functions 
- a necessary prerequisite to introduce the nonperturba¬ 
tive real-time functional RG. In Sec. m we discuss the 
properties of the functional flow and its application to 
calculate real-time response in equilibrium systems. We 
define the relativistic 0{N) vector model and associated 
scale-dependent effective action in Sec. m and discuss 
the role of different interaction vertices as well as con¬ 
sistent truncation schemes for the RG flow. The main 
part of this work consists of Secs. Ivl - IVIII where we in¬ 
vestigate different possible regimes for the dynamics, the 
properties of scaling solutions, and how our results com¬ 
pare quantitatively to known results from the e = 4 — d 
and large-iV expansion. We close in Sec. jVIIIj with our 
perspective on the quantum-to-classical transition in the 
low-temperature regime. 


The correlations and response of a system to external 
perturbations are fully captured by the time-ordered n- 
point functions: 


{f[(paAxi) ■ ■ ■ 

I 6^Z[J] 




Z[J] 


( 1 ) 


j=o 


where the Heisenberg-field operators (p°'(x) = 

given in the Wcomponent 
vector representation of the 0{N) symmetry group and 
the functional derivatives with respect to the real-valued 
external sources J°‘{x) are evaluated at spacetime 
points X = {t,x). For a system that is prepared in the 
groundstate or in thermal equilibrium, the generating 
functional in Eq. 0 can be defined as 


Z[J] = Tr|e-^^r 




■w]}. 


( 2 ) 


Here, /? = I/T defines the inverse temperature, H is the 
(interacting) Hamiltonian, and the operator T orders the 
fields on the real line. However, systematic perturbation 
theory based on 0 is ill-defined (see, e.g.. Refs. [43ll44] b 
To circumvent this point, we introduce the closed time 
contour C [33] [3S] , so that the generating functional can 
be written as 


Z[J] =Tr{e-^*Tc 





(3) 


The external sources are taken to have nonvanishing sup¬ 
port on parts of C in the complex time plane, upon 
which the operator Tc defines the time ordering of the 
fields. The contour integration in Eq. ^ can be ex¬ 
pressed via the parametrization C : [0,1] 9 s >->■ t{s) G C, 
• • • = ds F(s) f d'^x • • •, with t(0) = t(I) = U, 
and ti real. While the purpose of this section is to il¬ 
lustrate how to determine the real-time dynamic corre¬ 
lation functions and response, for now, we will impose 
no restrictions on the time path C. In the following, we 
examine different possible choices and highlight their im¬ 
portance and shortcomings. 

The imaginary-time formalism (ITF) [37IH3] provides 
a convenient means to calculate imaginary time-ordered 
correlation functions in the groundstate or in ther¬ 
mal equilibrium. In this framework, the helds live on 
an imaginary-time contour C = C-ip, where C-ifj = 
{t G C I Ret = ti, Imt G [—/3,0]} (see Fig. [^. That is, 
in the ITF the {d + l)-dimensional Euclidean theory ap¬ 
pears compactified in the time direction with perimeter 
fd = 1/T and the imaginary-time Green functions sat¬ 
isfy a periodicity constraint, the Kubo-Martin-Schwinger 
(KMS) condition [501 El]. Based on the generating func¬ 
tional ([^, this formalism allows us to determine thermo¬ 
dynamic properties of interacting many-body systems. 








[jL 
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FIG. 2. Imaginary-time contour C-ip in the complex time 
plane. The right branch defines the segment of the con¬ 
tour on which the sources are nonvanishing and is equipped 
with a imaginary-time ordering (indicated by the arrow). 
The left branch closes the contour and dehnes a periodicity 
constraint for the fields that live on the right contour (at 
initial and hnal imaginary time). In the corresponding Eu¬ 
clidean action this is expressed by the boundary condition 
(p“(t = ti) = v?“(f = ti — iP). The particular choice on which 
branch the fields should propagate is arbitrary. 




FIG. 3. Real-time contour C = C+ © C_ in the complex 
time plane. The upper and lower part of the contour define 
those segments on which the external sources are nonvan¬ 
ishing. The statistical operator po = is inserted 

at initial time ti, while the contour winds around at final 
time tf. The endpoints of the contour are completely ar¬ 
bitrary and the generating functional does not depend on 
the specific choice of ti ox tf. However, initial and final seg¬ 
ments define constraints for the helds that live on the upper 
and lower branches of the contour, e.g., in the real-time 
action on the closed time path, we impose the constraint 
^p\{t = tf) = (fi°L{t = tf) at final time. 


either by perturbative or nonperturbative means. The 
analytic continuation of imaginary-time correlation func¬ 
tions furthermore allows for the calculation of dynamic 
properties, e.g., the response to external perturbations 
and associated transport properties. However, this ana¬ 
lytic continuation, in practice, is an ill-posed numerical 
problem. In the absence of a method that produces ex¬ 
act results at imaginary time, additional assumptions on 
the analyticity properties and ultraviolet (UV) behav¬ 
ior of correlation functions are required (see, e.g.. Refs. 
[5^ [55] ). The validity of such assumptions are especially 
questionable in the vicinity of a thermodynamic singu¬ 
larity, e.g., at a second-order phase transition where the 
system displays nonanalytic behavior. Either way, the 
final result will be sensitive to the values of some ad hoc 
parameters that might not have an immediate physical 
significance. To determine the dynamic correlations or 
response, it is therefore desirable to employ a framework 
that allows for a direct calculation in real time without 
the need for analytic continuation. 

Real-time formalisms (RTF) [151 UU [531 [55] employ a 
different choice for the time contour C that appears in the 
definition of n-point functions 0 and their generating 
functional ([^. While in principle there is an infinite set 
of possible time paths, we consider a time-ordered con¬ 
tour consisting of two branches only, C = Cj^ ® C-, with 
t(0) = t(l) = ti and t{s) real for all s G [0,1] [5^ 155]. 
This contour, illustrated in Fig. [^ is defined with for¬ 
ward and backward time ordering on the branches C+ 
and C_, respectively. The closed time path (CTP) gen¬ 
erating functional 0 becomes 


where we have used the cyclicity of the trace and intro¬ 
duced external sources and J_ on the forward and 
backward time-ordered segments of the contour. We as¬ 
sume in the following that the sources have nonvanishing 
support at times ti < t < tf and fall off sufficiently fast 
at spatial infinity]^ This doubling of the degrees of free¬ 
dom allows us to employ the usual definition of the time 
ordering operator T for real times. Furthermore, writ¬ 
ing the generating functional in the form Q makes the 
conjugation symmetry [15] clear: 

Z[J+,J_Y = Z[J_,J+]. (5) 

Of course, in the absence of sources the physical content 
of Q is no different than that provided by evaluating the 
generating functional on the imaginary contour C-ip, and 
we have Z^"^^ = 

In the following, we aim at a nonperturbative evalu¬ 
ation of the generating functional using the functional 
RG in the RTF. Its construction relies on the path in¬ 
tegral representation of the generating functional on the 
CTP, the derivation of which [551I551I551I55] is based on 
the insertion of a complete set of states associated to the 
Heisenberg-held operator ip(x). We distinguish the helds 
on the two branches of the CTP and introduce two dis¬ 
tinct degrees of freedom, and (p-, on the forward and 
backward time-ordered segments, respectively: 

(p(t, x)l<fia,(x); tGCa) = x)\ipc,(x); teCa), ( 6 ) 

where a G {+,—}. At large time tf, at which we identify 
the two branches of the contour (cf. Fig.|^ the helds sat¬ 
isfy the property ip+{tf) = tp_(tf) and we will choose the 


Z[J+,J_] = Tr|T 




X e 




l\t 


, (4) 


^ The spacetime integration is defined as: fx ‘' = 

is understood that the limits 
ti —oo and tf oo are taken at the end. 
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value oi tf in the following such that t/ —>■ oo. Proceed¬ 
ing along these lines we first decompose the generating 


functional Q as a product of transition matrix elements 
between states that live on the CTP, i.e. 


Z[J+,J-]= n [d‘PaiU)][d(paitf)] {(p+-,U\e (7) 


where the subscripts J+ and J_ indicate that the cor¬ 
responding matrix elements should be evaluated in the 
presence of external sources. The functional measure is 
defined as J[d(pa{t)] = Af /na;eR‘^ d(pa{t,x) with some 
appropriate normalization constant Af and a proper reg¬ 
ularization of the infinite product. The matrix elements 
for which the external sources are nonvanishing can be 
expressed in terms of a functional integral over a phase 
factor which is fully determined by the dynamics and the 


coupling of the fields to external sources. On the 
a G {+,—}, segments they take the following form: 

Here, the functional integration runs over all field con¬ 
figurations with fixed values of the field both at initial 
‘Pa = Pa(ti) and final time pa = Pa(tf). Evaluating the 
product of all matrix elements in Q we obtain the path 
integral representation of the full generating functional: 


Z[J+,J_]= f 

aGf-h.-} 


(9) 


where we have introduced the CTP action 
S[p+,p-] = 5'b+] - S*[p_]+ 
and the full functional measure 

n = / n n idpait)]. 

f , 1 ^ f 1 la . 


( 10 ) 


( 11 ) 






The functional integral goes over all possible field config¬ 
urations. Any constraints on the fields have been ab¬ 
sorbed into the CTP action (10), in the form of the 


constraint-fluctuation (CF) functional J^[pj^,p_\. We 
provide the following formal definition in terms of the 
matrix elements at initial and final time 

T[p+,p-] = -i\nZ - i\n{p+-,ti\pD W-'.U) 

-i\n{p_-,tf\p+]tf), (12) 

where po = Z~^e~d^^ so that Trpu = 1. Of course, 
to give any meaning to this object, we need to consider 
regularizations of the contributions on the RHS. This 
is immediately clear if we examine, e.g., the final time 
constraint which can be expressed in the following way: 
{p-\tf\p+]tf) = 5[p-{tf) — p+(tf)]. Independent of the 
regularization we take note of the following properties: 
The CF-functional references only the field values at the 
boundaries ti and tf and, in equilibrium, it is homoge¬ 
neous in time. Also, while the first term on the RHS of 


Eq. (12) is field-independent and therefore does not con¬ 


tribute in the calculation of correlation functions the re¬ 
maining terms do. They eventually lead to nonvanishing 


mixing contributions and additional couplings between 
the p+ and p- fields in the CTP effective action. It is 
these terms that affect the correlations and response of 
the system and have to be accounted for to resolve the 
real-time properties of the theory in equilibrium. 

From Eq. (|^ we obtain the generating functional for 
connected correlation functions in the path integral for¬ 
malism. It is defined as 


1T[J+, J_] =-flnZ[J+, J_]. (13) 


Field expectation values are obtained by functional dif¬ 
ferentiation 


4>a{x) 


5W[J+,J_] 

SJa{x) 


as are the connected correlation functions 


Gap{x,y) 


S^W[J+,J.] 

5Ja(.x)5Jp{yY 


(14) 


(15) 


where Greek indices a, P G denote fields and ex¬ 

ternal sources on the distinct branches of the CTP. We 
construct the effective action in the RTF via the Legen¬ 
dre transform 


T[P+,cl>_] = W[J+,J_] 

- [ {J+A^)YI{X) - J-^a{x)(l)°L{x)} , (16) 

J X 

for which the sources are considered to be functionals of 
the fields, i.e., Jf. = J+[4>+A-] and J_ = J-[4>+A-\, 
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as obtained from inverting the relations (14). From the 
effective action we derive the equations of motion in the 
presence of nonvanishing sources 


5(j)a{x) 


- sgn(a) Ja(x). 


(17) 


When the external sources are set to zero, in the ground- 
state or in thermal equilibrium, the field configurations 
that extremize the effective action are homogeneous both 
in space and time. While in the construction of the ef¬ 
fective action we employed the and cj)- fields, a much 
more convenient parametrization of the physical degrees 
of freedom is instead given by the average and difference 
fields (retarded/advanced basis) 


4 ’=^{ 4 ’+ + 4 ’-), ^ = ( 18 ) 

and the corresponding external sources 

J=i(J+ + J_), J = J+-J_. (19) 


In this parametrization, the equilibrium field configura¬ 
tions are: I j=j=o = (l^-\j=j=o = 

4’\,j=j^o ~ define the retarded propagator in 

the following way 


equilibrium as a function of a given set of microscopic pa¬ 
rameters in the action. That is, we do not show how ther¬ 
mal correlations are built up in real time starting from 
an arbitrary initial state. This is an interesting prob¬ 
lem in its own right and there are other frameworks that 
address the issue of thermalization [57MT] . We assume 
that the averaging over the initial statistical operator has 
been taken into account to the effect of imposing rela¬ 
tions between different n-point functions dictated by the 
properties of equilibrium fluctuations. Typically, these 
relations are expressed in frequency space and this im¬ 
plies that the statistical operator po has been inserted 
in the infinite past ti —>■ — oo, so that the system satisfies 
time translation invariance. 


The derivation of the CTP action (10) makes it clear. 


that the calculation of correlation functions will not de¬ 
pend on the field-independent contribution —i In Z to the 
CF-functional. In fact, it can be ignored for the most 
part when dealing with a system in equilibrium and this 
is what we will do in the following. However, this proce¬ 
dure leads to some subtle issues in the calculation of the 
thermodynamic free energy, which we address here: The 
free energy density in the RTF is defined as 


F = -V - ^ - 

/3 /3VoWi 


InZ, 


( 22 ) 


G^ix,y) 


S'^W[J,J] 

6J{x)6J{y)’ 


( 20 ) 


which (in the absence of sources) is equivalent to the ex¬ 
pectation value of the anti-symmetric field commutator, 
i.e., G^{x,y)\^^^j^^ = i9{x°-y°){[p{x),<p{y)]). The sta¬ 
tistical correlation function is given by 


F{x,y) 


. 6^W[J,J] 
'sj{x)SJ{y)’ 


( 21 ) 


which, also in the absence of sources, can be 
written as the (connected) average anti-commutator 
F{x,y)\j^j^Q = ^{{(p{x),<p{y)}). While the derivation 
of the effective action relies on the Legendre-transform 
with respect to the external sources, in the end we are 
interested only in the case when they are set to zero. To 
avoid a proliferation of notation in the following, it will 
therefore be understood that all quantities are specified 
in the absence of sources unless stated otherwise. 

Equation ([^ essentially contains two different func¬ 
tional averages: The conditional average over the real¬ 
time evolution of the fields at times t > ti, as well as 
an averaging over the groundstate or statistical opera¬ 
tor pd at initial time ti. These two contributions are 
not distinguished within the ITF, where the functional 
average is performed over all (imaginary-time) fluctua¬ 
tions at once. Thus, in principle, the RTF allows us to 
study separately the effect of the initial conditions and 
the dynamics. Here, however, we ask specifically about 
the scale dependence of different dynamical processes in 


where Vold+i = dt f d'^x and 

Y = r[(/)= {ip),^ = 0] 

Vold+i 


(23) 


is the CTP effective potential. It depends on the value of 
the homogeneous field expectation value (ip) and we have 
pulled out the field-independent contribution ~ In Z from 
the CTP effective action to make it explicit. In equilib¬ 
rium, one finds the effective potential V to be zero in a 
calculation where all zero momentum IPI diagrams are 
summed over [331 [33]. Of course, this is simply a con¬ 
sequence of the conjugation symmetry ([^, which con¬ 
strains the free energy to be real-valued. However, it is 
important to realize that derivatives of the CTP effective 
action might be nonzero even in the absence of sources. A 
careful analysis shows that they match the contributions 
from the real part — ^ void+i ^ |b2j . Thus, the deriva¬ 
tives of the effective action carry the same information as 
In Z and the free energy can be reproduced also within 
the CTP formalism when the field-independent contri¬ 
bution to the effective action is not taken into account 
explicitly. 

The presence of conjugation symmetry ([^ constrains 
the form of possible parameters and couplings of the the¬ 
ory. For the CTP effective action it takes the form 

r[</),0]* = -F[0,-0]. (24) 

If we expand r[(^, (jj] in the fields (j) and 0 
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, Xm\ 2/1, • • ■ , yn)(l){xi) ■ ■ ■ (l>{Xm)Hyi) ’ ’ ’ ^(//n), 




r(’"’”^(a;i,...,a;m;2/i,...,2/n) = 


(5"^+"r[0, (^] 


(50(a;i) • • • S(j){xm)54>{yi) ■ ■ ■ 5(j>{yn )' 


(25) 

(26) 


where the summation over internal indices is implied, the 
symmetry (24) divides the contributions to the effective 
action into two distinct classes, characterized by either 
real or imaginary (amputated) IPI correlation functions 
T<^'^’^^Xi,...,Xm ; 2 / 1 , ■ • ■, 2/n) (which are related to the 


odd and even powers of (j) in ( [2^ , respectively). In fact, 
already from the form of the CF-functional ( |I^ , it is 
clear that imaginary couplings and parameters to the 
CTP action are necessarily present when initial fluctu¬ 
ations and boundary conditions are taken into account. 
They play an important role over a wide range of scales 
and ensure the consistency of the thermal flow by pro¬ 
viding relations between real and imaginary vertices of 
the Wilsonian effective action. The general aim of this 
paper is to demonstrate how these vertices might be gen¬ 
erated and to show how an effective classical description 
emerges in the RTF when all fluctuations have been ac¬ 
counted for. 


III. NONPERTURBATIVE FUNCTIONAL 
RENORMALIZATION GROUP IN THE 
REAL-TIME FORMALISM 


The nonperturbative RG employed in this work is 
based on the flow of generating functionals for correla¬ 
tion functions. In particular, for systems with spatial 
and time translation invariance, it is useful to consider 
the functional flow of the scale-dependent IPI effective 
action |20| . There exists a large number of reviews in 
the literature and we will not provide all the details here, 
see, e.g.. Refs. [6dH68] . However, we will highlight some 
of the properties that appear in the physical (real-time) 
representation. 

The construction of the functional flow equation re¬ 
lies on a modification of the original CTP action, i.e., 
0\ S'[(p, (p] + AfeS'[(p, 0\, where 


tp] = - / ip°‘{x)Rk,abix,y)p>^iy). (27) 

Jx,y 


The function Rk is taken to be nonnegative and its 
properties are chosen so that (27) defines a momentum- 
dependent mass term, i.e. 


Rk,abix°,y°;p) = Rk{p)S{x° - y°)Sab, (28) 

where ti < x°,y^ <tf. It depends on an additional scale 
parameter k which controls the effective mass of differ¬ 
ent modes. Thus, by an appropriate choice of 7?^, we 
may regulate the infrared (IR) divergences associated to 


r 


the massless modes in the vicinity of a second order (or 
weakly first order) phase transition (see, e.g.. Ref. |69)i. 
Further regulating contributions quadratic in the differ¬ 
ence field ip might be necessary in the presence of ad¬ 
ditional symmetries [SD], but this is not the case here. 
For our purposes, it is sufficient to consider a regula¬ 
tor coupling the p and p fields only. Such a regulator 
is compatible with the symmetries in the presence of a 
fixed reference frame - the rest frame in which we define 
the temperature of the system. 


The modification of the action by the insertion (27) 
endows the generating functional for connected correla¬ 
tion functions VF[J, J] with a scale dependence, which we 
highlight by an additional index k: 


Wk[J,J]=-i\nZk[J,J]. 


(29) 


Any n-point correlation function derived from Wk will 
inherit this scale dependence and by functional differen¬ 
tiation with respect to the external sources, we obtain 
the field expectation values 


4>k{x) = 


SWk[J,J] 
6J{x) ' 


^k{x) = 


5Wk[jJ] 

6J{x) 


(30) 


in the presence of nonvanishing sources J, J, and Rk- 
Similarly, the (connected) retarded propagator and sta¬ 
tistical correlation function read 


Gkix,y) = 


Fk{x,y) = -i 


S^WkjjJ] 
6J{x)SJ{y)' 
S^WkiJj] 


(31) 


(32) 


5J{x)6J{y) 

The advanced propagator is not independent and given 
by the relation G^{x,y) = [G^{y,x)] , while the anom¬ 
alous propagator Fk{x, y) = —i vanishes in the 
absence of sources and therefore does not yield any con¬ 
tribution to the RG flow. 

The scale-dependent generating functional of IPI cor¬ 
relation functions F k [4>, <! >] is defined by the (partial) Le¬ 
gendre transform of Eq. (29), with respect to the sources 


J and J, for fixed (i.e., fc-independent) values of the fields 
(j) and (j>: 

= Wk[J,J] - AkS[<l),^] 

- J I Ja(a;)0“(x)-h Ja(a:)^“(a;)| . (33) 

Here, J = J[4>, (p] and J = J[(j), p] are in fact fc-dependent 
and are determined by inverting relations (30) for fixed p 














and (j)- From Eq. (331 it is straightforward to derive the 


exact RG flow equation: Taking the scale derivative with 
respect to the logarithmic scale parameter s = ln(fc/A) 
(where A refers to some fixed reference scale, which will 
typically be identified with the UV cutoff), we obtain 




Ffc = i 


Tr<; ^Rfe(p)ReG^(a;,p) 


(34) 


where the integration runs over real frequencies and mo¬ 
menta, p = (uj,p), and the trace Tr{---} is evaluated 
over internal field indices. Both the scale derivative of 


the CTP effective action, on the LHS of Eq. (34), as well 


as the retarded propagator on the RHS are defined in 
the presence of arbitrary sources J and J and therefore 
depend on the fields. At this level the functional flow 


equation Eq. (34) corresponds to an infinite hierarchy of 


coupled differential equations. Taking functional deriva¬ 
tives with respect to the fields, we may derive the flow 
equations for n-point functions to arbitrary order. To 
solve this hierarchy however, we need to close the infi¬ 
nite set of coupled equations which defines a particular 
truncation of the scale-dependent CTP effective action 
Pfc. The so-obtained flow equations are then evaluated 
at the global minimum of the effective action, determined 
by (30) in the absence of sources, J = J = 0. Of course. 


any truncation of (34) amounts to an approximation of 


the full theory and restricts both the possible degrees of 
freedom and interactions that might become relevant in 
the low-energy effective theory. While sophisticated trun¬ 
cations have been developed to solve these flow equations 
(see, e.g.. Refs. [TOHZS]) typically the quality of a partic¬ 
ular truncation can be improved considerably if one has 
some understanding of the relevant operators in the low- 
energy regime. In the following section we will consider 
possible truncations for the 0{N) model that allow us 
to follow the flow of the scale-dependent effective poten¬ 
tial in the vicinity of the second order phase transition 
(at vanishing or nonvanishing temperature). This allows 
us to address the impact of the characteristic scales of 
the system determined by temperature and microscopic 
interactions and to resolve the dynamic scaling behavior. 

Let us finally comment on the form of the regulator 
function Rk{p)- So far it seems, that we might choose 
any function as long as it is nonnegative. However, if 
we are to make sure that Eq. (34) defines a functional 


flow between the microscopic action and CTP effective 
action, we need to impose additional constraints on the 
regulator function that enforce compatibility with 
these boundary conditions. In particular, the following 
properties should hold: limfe_>A Pfe = S[p,(p\ and 


^ By writing the flow equation in the frequency-momentum repre¬ 
sentation, we assume that the limits ti —>■ —oo and ty —> oo have 
been taken, thus restoring time translation invariance. Here and 
in the following we define the frequency-momentum integration 

Sp---= (2^)Vi ■ 


limfc_>o P/c<('] = <(>]■ Of course, when J = J = Q, 

the effective action should depend only on (f) = {(p) while 
(j) = 0. These conditions can be implemented by the 
following requirements: limfc_>A Rk ^ h? ^ oo, and 
limfc_,,o Rk = 0, respectively. Specifically, in the main 
part of this work, concerned with the critical dynamics 


at nonvanishing temperature (cf. Sec. VI), we employ the 


following frequency-independent regulator function 


Rk{p) = Z^{k‘^ - p^)e{k'^ - p2), 


(35) 


while at zero temperature we choose to work with a sim¬ 
ilar (d-|- l)-dimensional Euclidean regulator (cf. Sec. IVI). 

is a scale-dependent factor which appears in the def¬ 
inition of our truncation for the effective action Pj, (see 
Sec. IVB). Eq. (35) is also known as the Litim regula¬ 


tor and satisfies an optimization criterion that reduces 
the spurious scheme dependence from the truncated flow 
equations HZ]. Compared to other regulator functions it 
has the clear advantage that it allows us to derive ana¬ 
lytic expressions for the nonperturbative RG functions. 
However, a word of caution is in order at this point: Dif¬ 
ferent choices of regulators might affect the quality of 
the final results. To understand the possible effects that 
such a choice might entail, we compare the performance 


of different regulators in Sec. VH That is, we apply an 
expansion around the upper critical dimension dcr = 4 
for the classical phase transition and compare to known 
results from the e-expansion. A similar analysis is pro¬ 
vided for the large-fV expansion. This will give an indi¬ 
cation to what extent results obtained within the func¬ 
tional renormalization group are regulator dependent for 
a given truncation of the scale-dependent effective action. 

Before we proceed, let us remark that we will drop the 
index k from now on to avoid a proliferation of notation, 
unless it is necessary to highlight the scale dependence 
of parameters or couplings of the model. 


IV. RELATIVISTIC 0(A) MODEL 

In this section, we discuss the general structure of the 
CTP effective action and its truncations with the ob¬ 
jective to smoothly connect the coherent and incoherent 
dynamics in the UV and IR, respectively. In Sec. IV A 
the CTP action at T = 0 is discussed in detail, whose 
structure governs the microscopic coherent dynamics at 
arbitrary T. In Sec. |IVB[ our ansatz on the two-point 
propagators is explicitly constructed so that the fluctua¬ 
tion theorem holds exactly at the linear level. This fixes 
the scaling dimensions of the fields, which we provide 
in Sec. |IV C[ Finally, in Sec. |IVD| we discuss consistent 
truncations of the vertex expansion of the CTP effective 
action. These considerations allow us to address the ef¬ 
fectiveness as well as the limitations of the assumptions 
implicit in our ansatz. 
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A. Microscopic theory (T = 0) 


Here, we provide the CTP action for the relativis¬ 
tic 0{N) theory at zero temperature. We employ the 
((p+, (/?_) basis on the C+ © C_ contour (see Sec. 0 and 
start with the noninteracting theory S = Sq, for which 
we derive the propagators and discuss the implications 
of the CF-functional (cf. Eqs. (10) and ([HI))- This sec¬ 
tion serves essentially as an introduction to the struc¬ 
ture of the 0{N) theory in the real-time formalism (see 
also Ref. m which provides an extensive introduction 
to real-time methods in field theory). The propagators 


are defined and their properties outlined. While most of 
them generalize also beyond perturbation theory, care is 
taken to clearly highlight these properties and distinguish 
them from the assumptions valid only in the perturbative 
regime. At the example of the free theory, we argue that 
the effect of the boundary constraints and initial fluctu¬ 
ations can be fully captured by introducing an effective 
CTP action. The free action can be written in terms of 
contributions from the two distinct branches of the con¬ 
tour plus an additional contribution that takes the form 
of a mixing term ^ ip_: 




^ + 10+ — a;^(p) 


;)-l) 


-f0+(sgn(w) + 1) 

— — *0+ — w^(p)) 


( ^+{P)\ 


(36) 


where ui{p) = ± + inn? defines the characteristic fre¬ 

quency of the modes. Note, that the spectrum of charac¬ 
teristic frequencies and propagators of (36) are invariant 


under 0{N) symmetry (internal indices are not explicitly 
written out, unless necessary). This is by no means obvi¬ 
ous: In fact, a naive doubling of the degrees of freedom on 
the CTP would yield a theory which is invariant under in¬ 
dependent 0(A^)-rotations of the ip+ and p- fields. How¬ 
ever, such a symmetry is unphysical - the fields on the 
upper and lower branch of the CTP are not independent 
due to the presence of the boundary conditions that enter 
in the construction of the path integral, cf. Sec. We 
include the nondiagonal mixing contributions ~ 
in Eq. ( |36[ ) precisely to account for the nonvanishing cor¬ 
relation of the fields and to constrain the spectrum of the 
theory]^ Clearly, their effect is identical to imposing the 


® This situation is very much reminiscent of the situation that one 
encounters in the replica approach to spin glasses ED. As an 
example consider two identical copies of a classical lattice Ising- 
spin model { 5 “}, where the Greek index a labels the different 
replicas, while the Latin indices i,j label the sites of the discrete 
spatial lattice. The Hamiltonian, H = 2 j ~ 

£ is defined via the exchange couplings which are 

not necessarily nearest-neighbor and randomly alternate in sign. 
It includes an additional mixing term which introduces nonva¬ 
nishing correlations between replicas. Depending on the spe¬ 
cific choice for the probability distribution of the couplings, 
one might observe a large degeneracy of equilibrium states for 
such a spin system at sufficiently low temperatures. The cor¬ 
relations between different replicas define an order parameter 
low-temperature glassy phase m- 
The existence of one or more (coexisting) equilibrium states can 
be monitored by observing if is continuous or discontinu¬ 

ous when e changes sign m, and the presence of a discontinuity 
is related to the spontaneous breaking of the replica symme¬ 
try q°^^{0) = g^“(0) |82H84| . In the context of the CTP ac¬ 
tion the replica symmetry breaking is expressed on the level of 
the frequency-dependent propagators, i.e., 

e {+, 


I 

boundary constraint = tf) = ptit = tf) at final 

time tf. 

The presence of the infinitesimal contributions fur¬ 
thermore ensures the proper causality structure of the 
propagators and enforces a relation between the posi¬ 
tive and negative frequency modes. In the given field- 
representation the free propagators read: 


where 


G{p) = 


G++ip) = 


G++(p) G+_(p) 
G_+(p) G__(p) 


-1 


a;2 — LO^(p) + i0+ ’ 

G+_(p) = 2'Ki9{u))5{uj‘^ — u!^{p)), 
and the following symmetry relations hold: 

G++(p) = -[G__(p)]t, 
G+_(p) = G_+(-p). 


(37) 

(38) 

(39) 


(40) 

(41) 


The diagonal elements of the propagator (37) correspond 


to the usual time-ordered (and reverse time-ordered) 
Green functions. However, not all of the components 


are independent. Apart from the relations (40) and (41) 
the following algebraic identity holds 

G++(p) + G__(p) = G+_(p) + G_+(p), (42) 


as one may easily verify. Eqs. (40) - (|42[) reduce the num¬ 


ber of independent degrees of freedom and are valid also 
beyond perturbation theory. Given these relations one 


might argue that the propagator (37) can be expressed in 


terms of three independent real-valued components only. 
However, the presence of a fluctuation theorem [sniiHS] 
and the dispersion relation [86] further reduces the num¬ 
ber of independent degrees of freedom. Both of these re¬ 
lations are most easily derived in the retarded/advanced 
(RA) basis. Therefore, as we have done already in Secs.[H| 
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and III we will switch to this basis, where the properties 


of the action and the consequences of the iO'*" insertion 


become somewhat more transparent. In that basis, the 
free action takes the following form 




(43) 


Here, we see that the infinitesimal regulating contribu¬ 
tion appears as a positive-definite imaginary contribution 
~ to the action. Not only does this guarantee the con¬ 
vergence of the generating functional, but it also ensures 
that the dynamic correlations of the system are related 
to the response with respect to an external perturbation. 
In the RA-representation, the free propagator is given by 


G{p) = 


_ fiF{p) G^{p) 


\G^ip) 0 


(44) 


where 

G^{p) = 


-1 

a;2 — Lo'^{p) 


■ in sgn{uj) (5(a;^ — a;^(p)), (45) 


F{p) = n6{uj^ 


(46) 


define the free retarded propagator and statistical corre¬ 
lation function, respectively. By conjugation symmetry 
(24) the advanced propagator G^ is seen to satisfy: 


G^ip) = [Gi'(p)]^ 


(47) 


All components in ( |44[ ) are either purely real or imag¬ 
inary, i.e., = G^/^(—p) and [iF(p)]'l' = 

—iF{—p), in contrast to the propagators in the (p+, P-) 
basis. This is one of the advantages of expressing the ac¬ 
tion in the form (43) - it clearly emphasizes the relevant 


degrees of freedom. Furthermore, from the propagators 
(45) and (46) we may easily read off the following relation 


F{u},p) = sgn(a;) ImG^(a;,p) = 


P{^,P) 

2|a;| ’ 


(48) 


where in the last step, we have defined the spectral den¬ 
sity p{uj,p). In fact, this relation is exact - it holds for 
any theory that is prepared in the groundstate {T = 0). 
Eq. ( [4^ is a fluctuation theorem which relates the corre¬ 
lation for zero point fluctuations F{uj,p) to the spectral 
density (or alternatively the imaginary part of the re¬ 
sponse function). When the system is in its groundstate, 
all it can do is to absorb energy and without any cou¬ 
pling to a dissipative environment, the only possibility 
for non-zero fluctuations is to allow for virtual particle- 
antiparticle pair creation (and annihilation). This point 
of view can is supported by the following derivation which 
starts from the relation 


/ ia+\u}\f Faiuj,p) >0, 
J U) 


(49) 


where the correlation function is expressed in the mass 
eigenbasis, Fabip) = Fa{p)5ab- The components Fa{p), 
a = 1,... ,N, are real which follows from the reflection 
symmetry [E(p)]’f = F{—pY while positivity follows from 
the stability of the action|j Exploiting these properties, 
we may write Eq. (49) which is obviously satisfied for any 
real-valued parameter a. In particular, it holds true at 
the minimum (with respect to a) for which we obtain 

|y F;(a;,p)||y a;^Fa(a;,p)| > |y \u}\Fa{u},p) 

We use the following identities 

= f Fa{uj,p) 

J tJ 

^ J X 

^T^I{P)= f UJ^Faiw,p) 

J uj 


(50) 


(51) 


(52) 


where n(x°,x) = ^(p(x^,x} defines the conjugate mo¬ 
mentum to (fi. Together with the fluctuation theorem 
(48) we write 


y iwiFa(w,p)=^y Pa(uj,p )=y 


(53) 


In the last step of Eq. ( [53| we have employed the sum 
rule for the spectral density (see below), which follows by 
using the definition of the spectral function 

[ Paiuj,p) = -l f {[(Paix°,x),naiy°,0)])\^„^^„e-^P'^, 

J UJ ^ J X 

(54) 

and the canonical equal-time commutation relations for 
the field and conjugate momentum operators. This al¬ 
lows us to express Eq. (|50| in the well-known form: 


A^l{p)Anl{p) > 


1 


(55) 


® In the mass eigenbasis, where the propagators are diagonal, the 
statistical correlation function can be written as: 

Fa{p) = GR(p) " Vy- - . > 0, 

where F is the effective action (and corresponds to Sq in the free 
theory). 
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Finally, using the above expressions for the free prop¬ 
agators (45) and ([4^ one may easily check that the 


uncertainty relation is in fact saturated in the ground- 
state, where Aip^{p) = l/{2uJa{p))- Thus, for any finite 
momentum the groundstate features nonvanishing zero- 
point fluctuations, which shows that Eq. (48) is satisfied 
nontrivially. 

For completeness, we also provide the following fre¬ 
quency integrals: 




(56) 


that follow from the canonical equal-time commutation 


relations. In general, Eqs. ( 

53 

), ( 

55 

), and the normaliza- 

tion of the spectral function (56 

) 

lold true also in the 


interacting theory. This applies both to the bare and 
renormalized fields, for which the kinetic term has coef¬ 
ficient one in the effective action^ 

Apart from the fluctuation theorem, there is a disper¬ 
sion relation that relates the real and imaginary parts of 
response function. It is a direct consequence of analyt- 
icity of the retarded (advanced) propagator in the upper 
(lower) half of the complex plane 


^ReGi'(cc,p) = -P 


ImG^(a;',p) 


(57) 


where P denotes the Cauchy principal value. The reader 
may easily verify that this relation is satisfied by the free 
propagators, but it also holds for the exact propagators. 
Thus, it is clear that the propagator (44) of any system 


prepared in its ground state can be described in terms of 
a single real-valued function. In fact, this statement is 
also true at nonvanishing temperatures (see Sec. IVB|. 


This greatly simplifies the diagrams that one needs to 
evaluate in perturbative calculations and similarly, in the 
framework of the nonperturbative functional RG. 

Let us move on to the interacting case. The vertices 
in the CTP action S = Sq + Sx derive directly from 
the construction of the theory in the RTF illustrated in 
Sec. |n] For the Xip'^ theory, we obtain two vertices in the 
physical RA-representation, parametrized by the same 
coupling A, but different prefactors: 

Sx[ip,0\ =-^ j 

+ \p''{x)^a{x)(p^[x)^. (58) 

Their origin derives from the field degrees of freedom 
that propagate on the two distinct branches of the time 
contour C (cf. Sec. 0 which can be seen explicitly, by 


changing the field representation: 

Sx[p+,P-\=-j^ j {ipX{x)-ipt{x)] . (59) 

The additional minus sign originates from the reverse 
time-ordering along the C_ path. Thus, the two cou¬ 
pling terms in the RTF are simply inherited from the 
doubling of the held degrees of freedom on the C+ © C_ 
time integration contour. 

In the presence of interactions, both the mass parame¬ 
ter and the coupling of the theory are renormalized. This 
is due to the coherent interaction of bare particles with 
excitations from the vacuum and can be compensated 
by absorbing the associated UV divergences in appropri¬ 
ate counterterms. At nonvanishing temperature, how¬ 
ever, we expect an additional IR renormalization which 
is due to incoherent collisions with renormalized quasi¬ 
particles. These two processes lead to two different scales 
in the problem, which are characterized by the scattering 
length and mean free path of the system. The purpose 
of the following two sections is to illustrate the effect of 
these competing processes on the basis of an appropri¬ 
ate truncation of the CTP effective action. Eqs. (43) - 


(58) define our microscopic model that we will use as a 


starting point for the RG investigation. 


B. Low-energy effective theory 

In the framework of the nonperturbative functional 
RG, we employ an ansatz for the scale-dependent CTP 
effective action, which corresponds to a suitable trunca¬ 
tion of the gradient expansion and an expansion in field 
monomials (see, e.g., Refs. [87l|88]). Projecting the flow 
(34) onto those contributions included within our ansatz, 


we obtain a hnite set of RG equations for the correspond¬ 
ing parameters and couplings. This necessarily consti¬ 
tutes an approximation of the otherwise exact flow equa¬ 
tion. Here, we provide the specific form of the truncated 
CTP effective action in the 0{N) model and discuss the 
properties of the physical spectrum and dynamical be¬ 
havior as a result of the finite truncation. Fluctuations 
in the equilibrium state are reflected in the CTP effective 
action and it will therefore depend on temperature. We 
consider different temperature regimes and discuss how 
they are reflected within our ansatz. 

In the real-time formalism, the CTP effective action is 
most conveniently expressed in the frequency-momentum 
representation. We work in the RA-basis and assume 
that the effective action takes the following form 


^ However, note that certain truncations of the effective action (en¬ 
countered in the framework of the functional RG, cf. Sec. II V Bl 
might not retain the full spectrum of the theory, and in these 


cases the sum rule l |56[ l has to be modified accordingly. In partic¬ 
ular, this depends on the quality of the employed approximation 
for the frequency-dependence of the self-energy. 
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1 f ( 0 zW - + inoj/p\ fHp)\ 

2 Jp' ' yZ^^uj"^ — Z-^p"^ — i^luj/13 i(r2w//3) coth(/3a;/2) J\^{p)J 




(60) 


to second order in the gradient expansion. In contrast to 
the spatial momenta, an infinite series of frequency terms 
is taken into account, which is necessary to implement 
the fluctuation-dissipation theorem [50l[85] at nonvanish¬ 
ing temperature exactly at the level of the propagators 
(cf. @). The systematic construction of more elaborate 
truncations of the equilibrium effective action in the RTF 
significantly benefits from symmetry considerations |89j . 
In particular, one may show that fluctuation relations be¬ 
tween different (m,n)-point functions, = 4-, —rr-, 

\ I / ± ! d(p'^d<p'^ ‘ 

that characterize the properties of the equilibrium state, 
follow as Ward-Takahashi identities associated with time- 
reversal symmetry [89ll91j . 

The generalized potential lA = U[4>,4>) includes terms 
that appear to lowest order in the frequency and momen¬ 
tum expansion, while the wavefunction renormalization 
Z-^ and the renormalization factor Z'^ parametrize those 
terms that enter at quadratic order in momentum and 
frequency, respectively. The coefficient VL measures the 
strength of the imaginary contributions to the IPI two- 
point functions and All renormalization fac¬ 

tors Z^, Zll, and are assumed to be field-independent 
but depend on the RG scale parameter s. The consis¬ 
tency of Eq. ( [^ requires that > 0, /3 > 0, a con¬ 
strained imposed by the stability of the effective action, 
and Zy Z-^ > 00 is also scale dependent and 

defines the mass parameters and couplings of the theory. 
That is, the spectrum of the theory is given in terms of 
the N eigenvalues of the mass matrix squared 0 


m 


2 _ 
ab — 


— lim r 

p—>-0 


(Fi) 

ab 


(p) 


min 


jfu_ 

d4>°-dip^ 


(63) 


where the second functional derivative is evaluated at the 
global minimum of the effective potential: 


dV 


min 


limr(0’i)(p) =0. 

p—>-0 min 


(64) 


In equilibrium, the field configuration at the minimum of 
the effective potential is homogeneous and, without loss 


^ Note, that our definition of the renormalization constants is dif¬ 
ferent from the standard convention, i.e., Z-*- is defined as the 
inverse of the field renormalization defined in Refs. |92l I93| . 

^ Here, we define the two-point functions in the frequency- 
momentum representation with an appropriate factor for nor¬ 
malization: 


r(Fi)(p) = (r(M)(_p))t = 

Volti+1 

rlo.alfp) = (27r)^+i 5'^T[<t>,4>] 

Void_|_i Sti){-p)S4>{p) 


S^m4,] 

54>{-p)S(p{p) ’ 


(61) 

(62) 


r 


of generality, we assume that the field expectation value 
points in the 1-direction: 


= vS, 


al 5 


= 0 . 


(65) 


The propagators are most conveniently expressed in the 
basis, where the mass matrix is diagonal, i.e., = 

rn^Sab, and the two-point functions are diagonal. In par¬ 
ticular, we have 


-(1,1) 


(P) 


— Z-^uj'l{p) — i^ujjP 


^ab'i 


( 66 ) 


with the characteristic field-dependent frequencies 
Wa(p) = 3z^p'^ + a = 1,..., V, expressed in 

terms of the renormalized masses = rri^lZ^. The 
two-point function is exactly zero, while 


r 


( 0 . 2 ) 

ab 


(P) Lin = «(f^a;//3) coth(/3a;/2)<5ab, (67) 


at the minimum of the effective potential (65). Within 


our truncation there are no field-dependent contributions 
to r(0’2). This is enforced by two requirements: (1) 
that the linear frequency contribution should be field- 
independent and (2) that a fluctuation-dissipation theo¬ 
rem should hold. Before we comment on the properties 
of fluctuations and consistent truncations of the CTP 
effective action, we provide the general form of the prop¬ 
agators within our ansatz (60). In the basis where the 


mass-matrix is diagonal, the retarded propagator reads 


GabiP) = 


-1 


Zloj^ - Z^loUp) - R{p) + inuj/4 


Sab, (68) 


and the statistical correlation function takes the form 


Fab{p) = 


{p.uj/(3) coth(/3tLi/2) 


[ziuj^ - z^ioiip) - r{p)Y + {nuj/(3y 


(69) 

In contrast to the free propagators, they include the 
momentum-dependent function R{p) in the denomina¬ 
tor which provides a regularization of massless modes. 
The propagators corresponding to the IPI CTP effective 
action are obtained only in the limit when the RG scale 
parameter is removed k = 0 and all modes have been 
taken into account, i.e., R{p) = 0. 

With the given form of the scale-dependent propaga¬ 
tors it is easy to check that the following fluctuation- 
dissipation theorem (FDT) is satisfied: 

Fab{uj,p) = I ^n(w) -b 0 Pab{u},p), (70) 

where n{uj) = (e^“ —1)“^ is the thermal occupation num¬ 
ber, and the spectral density is given by 

Pab{u},p) = 2a;ImG^(,(a;,p). 


(71) 
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Let us briefly comment on the properties of Eq. (70) and 
its relation to the fluctuation theorem stated earlier in 
Sec. |IV A[ As far as our setup goes, we consider a closed 
system at a nonvanishing temperature, i.e., it is under¬ 
stood that different macroscopic subsets of the system 
are mutually in equilibrium. In a thermal state, we have 
the possibility that a given subsystem might fluctuate 
and exchange energy or particles with another subsys¬ 
tem. The total energy is conserved in this process and 
thus, strictly speaking, there is no dissipation. However, 
when degrees of freedom outside a given subset are inte¬ 
grated out, their combined effect might yield a dissipative 
coupling to a (classical) random potential in the resulting 
effective theory (see, e.g.. Ref. |94|)- In a similar spirit, 
the successive integrating out of modes that underlies 
the Wilsonian RG [SS] provides a low-energy effective 
description of the dynamics, where dissipation appears 
naturally through thermal fluctuations. Our aim is to 
show how the relativistic 0{N) vector model acquires a 
dissipative coupling and in what way it determines the 
dynamic universality class when the system is tuned to a 
second-order phase transition. To understand the transi¬ 
tion from unitary to dissipative dynamics Eq. (60) should 
interpolate between zero-point fluctuations at the small¬ 
est scales and macroscopic thermal fluctuations that ap¬ 
pear only through the RG averaging procedure. This 
is achieved by taking into account the exact form of the 
fluctuation theorem ( [70| ) in the framework of the CTP ef¬ 
fective action. Indeed, taking the limit /3 = I/T —>■ oo, we 
see that the FDT reduces to the fluctuation theorem pro¬ 
vided in Sec. |IV A| and by construction, the linear FDT 
is implemented exactly within our ansatz. Similar fluc¬ 
tuation theorems also hold for higher n-point functions 
characterizing the nonlinear response of the system |96j . 
We will come back to the nonlinear version of the FDT 
in Sec. IV D[ when we construct consistent truncations of 
the functional RG. 

Equation (601 constitutes the simplest possible ansatz 
that allows us to understand the competition between 
coherent propagation of particles and dissipation. It is 
important to understand the limitations of such an ap¬ 
proximation. For that purpose let us take a closer look 
at the limits of the effective GTP action and the asso¬ 
ciated dynamic behavior. At T = 0 the GTP effective 
action (60) maps onto the free microscopic action (43), 
up to wavefunction renormalization factors and 
that satisfy = Z-^ = I in the noninteracting case. In 
the presence of interactions - when all fluctuations have 
been integrated out - the phase structure of the model is 
captured fully by the parametrization of the effective po¬ 
tential, while the dynamics is governed by N propagating 
modes 


cc = ±{Z^/Z\\) ^ ujaip) + ^ (D/Z«)/3-' +0 (r") • (72) 

In the low-temperature regime the characteristic frequen¬ 
cies are renormalized, i.e., to leading order: uiR^a(p) = 

±(z-L/zi)® Wa(p) + 0(/3 ^). A nonvanishing imaginary 


part corresponding to a finite dissipative width appears 
at subleading order in the low-temperature expansion. If 
D/Zll is sufficiently small and Z-^, ~ I, one may speak 

of well-defined single-particle excitations (quasiparticles) 
with a nonvanishing mass and finite lifetime. 

This is in contrast to the high-temperature limit, where 
we are concerned with the dynamics of collective modes 

uj = i{Z^/Vl)ujI{p)P + 0{P^). (73) 


Their renormalized characteristic frequencies are given 
by WR^a{p) = (^~‘~/D) a;a(p)/Ip^ We first consider the 
hydrodynamic regime |p|C ^ lin the disordered phase. 


where the correlation length ^ 


is finite. We ob¬ 
serve that the characteristic frequencies do not vanish 
in this regime - any long-wavelength excitation will re¬ 
lax locally to the thermal state, with the typical relax¬ 
ation time T = uj~l^ ~ Thus, the hydrodynamics 
in the disordered phase is described completely in terms 
of non-propagating relaxational modes. As we lower the 
temperature to its critical value, the correlation length 
^ diverges and the order parameter experiences critical 
slowing down r ^ (for p = 0), where z is the dynamic 
critical exponent. The value of this exponent depends 
critically on the presence of slow (massless) modes and 
on whether they couple to the order parameter HHIU]. 
While the presence of such massless modes is ruled out 
in the disordered phase (assuming that there are no con¬ 
served quantities that couple to the fields), this is not 
so in the symmetry-broken phase, where N — 1 Nambu- 
Goldstone modes govern the low-energy dynamics. It is 
these modes that will - in principle - couple to the order 
parameter and determine the dynamic universality class 
of the phase transition [nma. 

In effective theories for critical dynamics, this mode¬ 
coupling appears as a phenomenological input from hy¬ 
drodynamics [98) . However, within a first-principles ap¬ 
proach, that starts from the microscopic action and at¬ 
tempts to connect to the low-energy dynamics, it is clear 
that they must follow from general symmetry consider¬ 
ations, that is, the fulfillment of conservation laws and 
Ward identities [99l 1100) . Within the framework of the 
functional RG these identities are satisfied only with the 
accuracy of the employed truncation [lOIL 1102) and in 
that sense, our ansatz Eq. (60) defines a nonconserving 
approximation to the low-energy dynamics. 

We therefore expect that within our approximation the 
low-energy dynamics in the vicinity of the continuous 
phase transition should follow that of Model A[iii[im] 
which is defined in terms of a nonconserved order pa¬ 
rameter in the presence of thermal noise. This makes it 
clear that our truncation will not allow us to distinguish 
different types of relaxational behavior (as required, e.g.. 


Note, that another set of finite-frequency modes decouples in 
the limit of small /3 and does not contribute to the low-energy 
dynamics. 
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to identify the proper dynamic universality class for the 
given microscopic model). Nevertheless, Eq. (60) serves 
as a simple, well-controlled starting point to explore the 
RG flow in the RTF and to identify how dissipative dy¬ 
namics emerges in the effective description at large scales. 
We leave a systematic study of different truncations and 
their associated possibilities for the dynamic critical be¬ 
havior for future work. 


C. Scaling dimensions and dynamic scaling 
relations 


It is instructive to consider the dimensions of the fields 
and parameters that appear in the microscopic theory 
and their scaling behavior when fluctuations are taken 
into account. The scaling dimension = [$] of the 
(possibly composite) scaling field $ is defined in terms 
of its behavior under scaling transformations x —>■ e~’^x 
and t —>■ 


$(a;) = e-^x), (74) 


where s = Infc/A is the scaling parameter (as defined in 
Sec. Ill). $ could correspond to the (p+ and ip- fields (or 
their counterparts in the RA-basis) in the CTP action 
(36), or the corresponding effective degrees of freedom 
that appear in our ansatz for the scale-dependent CTP 
effective action, Eq. (601. In general, each of these fields 
might have different associated scaling dimensions. 


The scaling dimensions of the spatial and time deriva¬ 
tives follow from their properties under scaling transfor¬ 
mations: 


[0] = 1, [5o] = 2. (75) 

In the presence of spacetime symmetries the scaling di¬ 
mensions are not independent and the dynamic exponent 
z can be expressed in terms of the scaling dimension of 
the spatial derivative operator. In particular, at T = 0, 
Lorentz symmetry enforces z = 1 exactly for the rela¬ 
tivistic 0{N) vector model. However, to define a nonva¬ 
nishing temperature, we need to specify a fixed reference 
frame. This breaks the Lorentz symmetry and for a ther¬ 
mal ensemble we find z ^ 1. 

To determine the dynamic scaling exponent z, or the 
scaling dimension A$ of the scaling field $, is a difficult 
problem in general: It requires to solve for the spectrum 
of the dilatation operator, i.e. 

(A$ -|- ztdo + X ■ d) $(<, x) = 0. (76) 

In the absence of interactions and at zero temperature 
however, this is an easy exercise: Assuming for the mo¬ 
ment that the fields (f>+ and (jj- (or equivalently and 


(fi-) are scaling, we obtain their scaling dimensions 


A+=A_ = -(74-2), 


(77) 


where D = d + z, and z = 1 at zero temperature. The 
equality A_(_ = A_ follows from the constraint imposed 
by the fluctuation relation at T = 0 and does not apply 
for nonvanishing T (see Eqs. ( [m| - (86) below). In fact, 
for T 7 ^ 0, the field s (j)+ and (j)- do not correspond to scal¬ 
ing fields (see Sec. IV C 2 below) and the corresponding 


scaling dimensions A± are not well-defined. 

In the presence of interactions, the fields are renor¬ 
malized and the canonical dimensions are modified. 
This holds true even at zero temperature, where the dif¬ 
ference in the scaling dimensions is captured by a single 
multiplicative factor 


\l Z^(t>a = 4>R,a, (78) 

a G {-)-, —} and define the renormalized degrees of 
freedom corresponding to the bare fields (j)a (in the limit 
fc —)■ A, -G 1). Although the renormalization factors 

appear as a simple rescaling of the fields, it should be kept 
in mind that they are in fact associated to the renormal¬ 
ization of composite operators, i.e., ^ Z^{d(f)a) ■ In the 
RA-representation, we use the following definition of Z^ 

zi{d4>+f - z^{d4>_f = z^d^^ ■ d(t>a, (79) 
while Zll is defined as 

zl{do^+Y - zi(5o</>-)" = Z^^dordo^a- (80) 


1. Quantum scaling regime (T = 0) 

At T = 0, we have = Z-*- = Z by Lorentz symme¬ 
try. Here, our interest lies in the scaling behavior at the 
continuous quantum phase transition where the fields ex¬ 
hibit scaling. That is, in the IR limit 0 < fc/A <C I, when 
all fluctuations have been taken into account and the mi¬ 
croscopic parameters have been tuned appropriately, the 
wavefunction renormalization satisfies Z-*- ~ . In 

this quantum critical regime we may define the anoma¬ 
lous scaling exponent 

77^ = -^lnZ^, (81) 

and it is clear that Z^ ^ with yyll = rj^ = rj. From 

the multiplicative renormalization of the fields we find 


The microscopic fields (pot and effective fields (l)oL,k (in effec¬ 
tive theory at scale k) have the same scaling dimensions in a free 
theory, i.e., 9 ?^ = 
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that the scaling dimension of the renormalized composite 
operator d<pfi^a ■ is given by 

A + A = D-2 + ?7^. (82) 

This result follows from an appropriate rescaling of the 
fields 4>r and 4>r^ from which we deduce the scaling di¬ 
mensions: 

A = A = ^(Zl-2 + ry^). (83) 

As in the free massless theory they are degenerate, which 
follows simply from the properties of the propagators un¬ 
der scaling transformations 

G^(w,p) = e-"(^+^-^)G^(e""a;, e®p), (84) 

F{uj,p) = e*p), (85) 

and the presence of the T = 0 fluctuation theorem 

F{uj) = sgn(a;)ImG^(a;). ( 86 ) 


2. Classical scaling regime (T ^ 0) 


At T ^ 0 the characteristic fluctuations in the vicinity 
of the classical phase transition (CPT) show a different 
behavior that is dictated by the scaling form of the FDT. 
In particular, we observe that 


J^(a;) = -ImG^(a;), 

LU 


(87) 


in the scaling region |a;| ^ T. We might encounter such 
a scaling regime if T is much larger than the cutoff scale, 
i.e., 0 < |tLi| < A^ ^ T. Using the scaling form of the 


dimensions for the fie. 


propagators Eqs. (84) and (85), we derive the scaling 


ds in the RA-representation: 
A = A — 0 , 


( 88 ) 


Thus, we observe that the presence of a nonvanishing 
temperature lifts the degeneracy in the scaling spectrum. 
This is equivalent to the statement that the retarded and 
statistical propagators acquire different scaling proper¬ 
ties at the classical phase transition (and the difference 
is quantified in terms of the dynamic critical exponent 
z). In contrast, at the QCP, both the correlations and 
the response are described by the same scaling behavior 


(A = A). Although, Eq. ( 88 ) indicates that the fields 
must scale differently, it does not provide the explicit 
form for the scaling dimensions at nonvanishing temper¬ 
ature. In Sec. VI we show that relation ( 88 ) is solved 
by 

^ 2 + ? 7 '^), A = i(d-2-b77^)-bz. (89) 


characterized by a continuous transition in the scaling 
dimensions of the fields, parametrized by the effective 
dimension 

A = ^ (de - 2-b??-^(de)) , (90) 

i.e., from d^ = d -I- 1 (w T) to dg = d (w <C T). In 
the literature this is often referred to as dimensional re¬ 
duction. However, from our discussion it should be clear 
that this statement applies only to the scaling proper¬ 
ties of the theory - the real dimensionality of the system 
remains unchanged. 


3. Dynamic scaling relations 


Let us consider the dynamic properties at the continu¬ 
ous phase transition: In general two parameters, jZ^ 
and U/^i, are necessary to specify the dynamics of our 
model (cf. Eqs. (72) and (IZH)) , while the wavefunction 
renormalization Z^ contributes to the renormalized mass 
spectrum. At zero temperature the finite renormalization 
of the theory requires only a single independent parame¬ 
ter Z^ (while Z^ jZ^ = 1 and H = O"*') and there is only 
a single anomalous scaling exponent, i.e., r] = = vK 

that features in the scaling spectrum. In contrast, at 
T 7 ^ 0 we have up to three anomalous exponents that 
contribute to the dynamic scaling in the vicinity of the 
classical phase transition. They are given by: 


as 


OS 


(91) 


With these we obtain the dynamic critical exponent z 
from the scaling behavior of the retarded propagator (84) 
and the statistical correlation function (85). If Z ^, ZJ ^ 


and U 7 ^ 0 , then each of the corresponding anomalous di¬ 
mensions will contribute in the scaling regime. Using the 
scaling assumption, we may derive the following relations 
within our truncation 


Z* > 0 : -ryll -b 2z - 2 -b 77 ^ = 0, (92) 

U > 0 : -T]^ -b z - 2 -b = 0. (93) 

If both Zll and U are nonvanishing, we arrive at the fol¬ 
lowing scaling relation: 

2 -b 27y^ - ?y^ - ry* = 0. (94) 


Thus, in general, it is possible to characterize the dy¬ 
namic scaling in terms of two independent anomalous 
exponents only, which we choose to be ?y'‘ and rp". The 
presence of additional symmetries may further reduce 
the number of independent exponents that character¬ 
ize the dynamic correlations. In particular, in the pres¬ 
ence of Lorentz symmetry (U = 0'^) we observe that: 


To summarize, regarding the dynamic critical behav¬ 
ior, we may distinguish between two scaling regimes: 


We may therefore conclude that there is a crossover be¬ 
tween the quantum and classical scaling regimes that is 


r 1 , T = o, 

\ 2 - Ty-L -b 7y^ , T 7 ^ 0 . 


( 95 ) 
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At zero temperature, in the absence of a dissipative cou¬ 
pling n, the system is invariant under Lorentz transfor¬ 
mations, which implies: = r]-^ = r], and we find that 

the dynamic critical exponent z = 1 , independent of the 
local interactions (provided that they do not break the 
Lorentz symmetry). However, the presence of a nonvan¬ 
ishing temperature leads to an additional relevant param¬ 
eter, which leads to a dynamic scaling exponent z > 1 . 


D. Local interaction approximation 


Our ansatz for the scale-dependent CTP effective ac¬ 
tion relies on a truncation of an expansion in gradients, 
as well as on a truncation of an infinite series of vertices. 
Here, we address the properties of the finite vertex ex¬ 
pansion, which rests on two assumptions: 

(I) We assume that the generalized potential U = 
(j>) takes the form of a formal series expansion 
that is local both in space and time: 

BV 

W = + E (96) 

^ n>2 


where V = V{cj)) corresponds to the effective poten¬ 
tial and the coefficients Taia 2 ...a„ = •?aia 2 ...a„(?!>) 
define the fluctuation amplitudes. 


(H) We assume that a truncation of the series (96) at 
4th order in the fields provides a reasonable ap¬ 
proximation of our theory in the low-energy regime. 
In this case, the only nonvanishing amplitudes 


Tn. 


that contribute in the generalized poten¬ 


tial are given by: 


^ aia2 — ^0102^2,0 “t ^ai^a2^2,2: (97) 

^a\a2a2 ^(aia 2 ^ 03 )-^3,1; (98) 

*7~aia2a3a4 ^(aia2^a3a4)774,0- (99) 


Here, the brackets denote a complete 
symmetrization in the indices, e.g., 

^{axa2^a2) — 3 (<^0102 <('03 + ^0403'^^’02 + Sa^a^'Pai)- 


Typically one expects that such a finite expansion is 
valid to a good approximation in the vicinity of a continu¬ 
ous phase transition, where higher order contributions in 
the fields are irrelevant in the RG sense. In the real-time 
formalism however, this might not be sufficient to resolve 
the appropriate low-energy dynamics. The vertices might 
carry a frequency- or momentum-dependence and this in¬ 
formation needs to be taken into account to address the 
question of the relevance or irrelevance of hydrodynamic 
modes in the low-energy limit. Indeed, using a local se¬ 
ries expansion (96) in our ansatz for the CTP effective 
action we may only determine the dynamics of the or¬ 
der parameter, without any coupling to long-wavelength 
modes. Thus, within our truncation, we expect that the 


IR dynamics is purely dissipative and the dynamic uni¬ 
versality class of the theory is given by that of Model A 
m, in line with our arguments at the end of Sec. 

The presence of a fluctuation theorem in the ground- 
state or thermal equilibrium constrains the form of the 
generalized potential. One finds that in general the 
vertices are not independent (see, e.g.. Ref. [35] where 
the fluctuation relations are given on the level of the 
3- and 4-point functions). These nonlinear fluctuation 
relations take a simple form if one uses a local vertex 
expansion. Together with the conjugation symmetry 
U{4>, ^)* = —(^), we find at 4th order in the fields: 

^aia 2 = •7'aia2a3a4 = 0. However, we point out that 


IV B 


we obtain a nonvanishing contribution to 

the CTP effective action from the zero-frequency limit 
of the quadratic term (jp' in Eq. (60). The simplifica¬ 
tion on the level of vertices is essentially a consequence 
of our assumption of locality and the requirement of con¬ 
sistency that we impose on our ansatz. That is, in our 
construction of the IPI CTP effective action, we have 
chosen to incorporate the properties of the equilibrium 
state exactly by imposing fluctuation relations between 
(to -|- n)-point functions Any consistent trunca¬ 

tion of the full equilibrium CTP effective action should 
take into account these relations [M] and this is what we 
have done here. 

In the following it will be useful to consider an alter¬ 
native representation of the fields. In particular, we in¬ 
troduce the following invariants under 0{N) transforma¬ 
tions 


cr 2 = CT3 = (100) 

in terms of which U = U{ai,( 72 ,( 73 ). This choice of 
parametrization makes the 0(V)-symmetry of the the¬ 
ory manifest and proves to be especially convenient in the 
derivation of the RG equations. Using the constraints 
that follow from conjugation symmetry and nonlinear 
fluctuation relations, we obtain the following form of the 
generalized potential: 

U = nPa2 + \i,2 (cti - v'^/2) 02 + A2,30-2cr3, (101) 

which we expand around the scale-dependent minimum 
(Tilmin = In the disordered phase (symmetric 

regime), wp > 0 and = 0, while in the ordered phase 
(symmetry broken regime), rrp = 0 and > 0. Since 
the functional RG provides a regularization of IR diver¬ 
gences associated to massless modes, we can follow the 
RG flow of the theory, through the phase transition, from 
the symmetry broken into the symmetric phase. Thus, 
we will encounter both scenarios when we solve for the 
RG flow equations. 

Note that all parameters and couplings of our model 
are real-valued. A generic ansatz for the generalized 
potential might also include imaginary couplings, e.g., 
*Ai_ 3 cricr 3 (see our discussion at the end of Sec. HU). The 
reason that they are not present here, is a consequence 
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FIG. 4. = 1 scalar model {D — 3, T = 0): As we tune the (bare) quartic coupling Aa = Aa/A at a fixed value of the squared 

field expectation value v% = v\/K = 0.078, we observe a phase transition, when the renormalized mass m|j vanishes. We define 
m\ = ^rVr in the symmetry broken regime and = kv'^, Xn = kX. This is shown in the left panel where follow the RG flow 
of m|j as the bare coupling Aa is tuned through its critical value AA.cr- In the middle panel we display the scale dependence 
of the renormalized coupling. At the quantum critical point the correlation length ^ diverges and the theory exhibits 

scaling. This scaling at the QCP is governed by the IR stable Wilson-Fisher fixed point, for which we determine the scaling 
exponents rj ~ 0.055 and z = 1. The scale dependence of the anomalous dimension is shown in the right panel. If Aa = Aa.cf 
the anomalous dimension assumes a constant value in the scaling regime (fe —> 0). 


of the fluctuation relations. Furthermore, not all cou- and couplings Xm,n is exactly canceled. To highlight 

the fact that we have rescal ed th e field s in such a way to 
compensate for the scaling (|102|) and (103) in the quan- 


plings in Eq. (1011 are independent. We may establish 
a relation between Ai ^2 and A 2 , 3 . We demonstrate this 
relation explicitly at the example of the flow equations 
(cf. Sec.©. 


V. QUANTUM REGIME (T = 0) 

The general form of the flow equations within our trun¬ 
cation is given in Appendix [A| Here, we provide the RG 
flow equations specifically in the quantum regime. This 
section summarizes well-known results [501 110411105j and 
serves to provide a complete picture in the context of the 


quantum-to-classical transition (cf. Secs. VI and VIII). 


The flow equations are derived by using a Euclidean reg¬ 
ulator function after a Wick rotation to imaginary times 
and frequencies. We employ the (d-l-l)-dimensional Litim 
regulator function m- To solve the flow equations it 
is useful to apply our knowledge of the possible scaling 
solutions (cf. Sec. IIV C]) : We define dimensionless renor¬ 
malized parameters and couplings in terms of which it is 
convenient to identify possible fixed point (EP) solutions 
of the RG flow. Let us consider the scaling dimensions 
for the parameters and couplings that enter in our model 
at r = 0: 


turn critical regime, we denote the corresponding param¬ 
eters and couplings with a bar. We also define Cl = fl, 
Zi = Zll, and Z-L = Z-L even though these parameters 
have not been rescaled. 

In the absence of a temperature scale Lorentz symme¬ 
try requires z = 1 and the RG equations simplify consid¬ 
erably. In particular, we find that no dissipative coupling 
is generated by the RG flow: 


OS 


(106) 


Furthermore, no imaginary couplings are generated and 
therefore they do not enter in the low-energy effective 
theory. This holds true independent of the choice of 
our truncation. In the following we use the relation 
Zll = Z-^ = Z and define the anomalous dimension 
rj = —^\nZ also outside the critical region. Thus, rj 
will in general be scale-dependent, i.e., rj = rik and takes 
on its critical value only in the scaling regime. 

The RG flow equations for the squared field expecta- 


[v^]=D-2 + 7^^, 

(102) 


[\m,n\ =4-£)-2?7-^. 

(103) 

dv^ 

We may define the following quantities 


ds 

-2 ^ 

(104) 


\ — ( 

^7n,n — ^ 

(105) 

dXi^2 

that are scale-independent at the critical point. 

. This is 

ds 

immediately clear, since we have Z^ ^ in 

the seal- 


ing regime, so that the scaling behavior of the parameters 



)v^ 


Sd{v) 


(H-u2Ai,2)' 


-b 1 


(107) 


+ Sd{v)>H,2 


(1 -b u2Ai,2)^ 


-biV-lL (108) 
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FIG. 5. = 1 scalar model (T = 0): Gaussian FP and 

Wilson-Fisher FP as a function of dimension D = d + 1. Both 
fixed points merge at the upper critical dimension Dcr = 4. 


^ = (i^- 4 + 277)A2,3 

+ '5z>(r?)Ai,2A2.3 |-^—^ - l| .(109) 

Here, the factor 


Sd{v) 


1 

2D-i7rfr((H + 2)/2) 



V \ 

D + 2)^ 


( 110 ) 


originates from the spherical integration in momentum 
space for the contributing diagrams. Observe that the 
quartic couplings satisfy the following relation: 


T 9A2,3 t dXi^2 

- — ^2,3^:^-, 

OS OS 


( 111 ) 


which clearly shows that they are not independent. From 
the microscopic action (431 we see that: Ai ,2 = 4 A 2,3 = 
A/3. By an appropriate rescaling of A and by virtue of 
Eq. ( |111 ) this relation between the couplings holds true 
for all RG trajectories in the quantum regime (T = 0). 

To close the system of RG equations, we need to pro¬ 
vide the anomalous dimension: 


rj = SoiO) 


71212 

V Ai 2 


(1 -f u 2 Ai, 2 )^ 


( 112 ) 


We find two possible fixed points (imposing the require¬ 
ment that the squared field expectation value is posi¬ 
tive, and both couplings Ai _2 and A 2,3 are nonnegative): 

(I) Gaussian fixed point 

^ 1,2 = ^ 2^3 = 0 , ( 113 ) 


and 77 = 0 . 

(II) Wilson-Fisher fixed point (2 < Z? < 4): 

Ai^2 = 4A2,3 ~ "3 ^ (114) 

and rj determined by Eq. ( |112[ ) where Ai 2 and 
are set to their respective fixed point values. 



FIG. 6. = 1 scalar model (D = 3, T = 0): The crossover 

scale fc* that separates the canonical from the nontrivial scal¬ 
ing region is set by the quartic coupling at the cutoff scale 
A. There is a one-to-one correspondence between the scale k, 
(upper panel) to points that lie on the curve connecting the 
Gaussian FP with the Wilson-Fisher FP (lower panel). As we 
move along the phase boundary (indicated by the arrow), sep¬ 
arating the spontaneous symmetry broken phase (SSB) from 
the symmetric phase (SYM), towards the Wilson-Fisher FP 
the crossover scale is shifted to the left (upper panel). Here, 
the maximum correlation length ^ = mff takes a finite value 
since both the bare values of the coupling Aa and field expec¬ 
tation value vf have not been tuned perfectly to the phase 
separation line. 


The behavior of the RG flow can be illustrated by con¬ 
sidering the D = 3 scalar theory {N = I) as an example, 
cf. Fig. In this case both fixed point solutions are 
present and the theory admits nontrivial IR scaling in 
the Ising universality class. We obtain 77 ~ 0.055 at the 
phase transition. This result should be compared with 
77 = 0.03639(15) from the high-temperature expansion 
|106j which seems to be the most precise determination 
to date (see, e.g.. Ref. [53] for a compilation of results and 
corresponding references). The discrepancy is not very 
surprising, since we have chosen a very crude truncation. 
Typically higher orders in the gradient expansion, e.g., 
to 4th order [75117511551 [TUT] provide a better numerical 
estimate of the scaling exponents in the context of the 
functional RG. In iA = 4 dimensions the only fixed point 
is the Gaussian one and the static and dynamic scaling 
properties are characterized by the respective mean-field 
values: z = 1 and 77 = 0. The trajectories of both fixed 
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points in the space of parameters and couplings, as we 
continuously vary the dimension D = d + 1, is shown in 
Fig. 0 

The RG flow in the 0{N) model features a character¬ 
istic scale at which fluctuations become important. This 
is shown in Fig. at the example of the D = 3 scalar 
model (TV = 1): Starting from the microscopic action, 
the anomalous scaling exponent crosses over from the 
canonical scaling behavior with 77 = 0 in the vicinity of 
the Gaussian FP, to the Wilson-Fisher FP with i] ^ 0 
at some characteristic scale k^. This crossover scale de¬ 
pends only on the value of the dimensionless quartic cou¬ 
pling Aa = AaA^“^ and is to a good approximation in¬ 
dependent of v\. When the coupling Aa is small, e.g., 
for £) = 4 — e close to the upper critical dimension, the 
Aa dependence of A* is seen to originate from the per¬ 
turbative one-loop corrections. That is, in an RG treat¬ 
ment these perturbative contributions become large when 
the RG scale parameter k approaches fc*. In Fig. we 
see that this perturbative argument remains valid even if 
e = 4 — H is of order one. 


VI. CLASSICAL REGIME (T / 0) 


In the presence of a nonvanishing temperature we ob¬ 
tain a different set of flow equations. Here, we ex¬ 
amine the scenario where the temperature is compa¬ 
rable (or much larger) than the cutoff scale A, i.e., 
(3A^ = I/Ta < 1, and the thermal occupation of modes 
n{uj) ~ Ta{uj/A^)~^ ^ I, is strongly enhanced. In this 


regime, the FDT takes the form: 


F{uj,p) = 


Ta 

ujA-‘ 


Im G^(a;,p). 


(115) 


As in the zero temperature case, it is useful to ex¬ 
press the flow equations in terms of the dimensionless 
renormalized couplings and parameters, which allows us 
to easily identify scale-invariant behavior. However, in 
contrast to the scaling behavior given in Eqs. (102 1 and 


(1031, there is an additional scale dependence in the 


classical regime which is introduced by the temperature 
T ^ A^. By an appropriate rescaling of the fields (cf. 


Sec. IV GI we find that the scaling dimensions of the pa¬ 


rameters and couplings are modified. In particular, the 
rescaled field expectation value squared and couplings 
take the form: 


d^=v\k/Ay, (116) 

= A„,„(fc/A)(™+"-4)A (117) 


This rescaling is sufficient to eliminate the explicit k- 
dependence in the flow equations. Apart from the pa¬ 
rameters and couplings in Eqs. (116 1 and (117), we de¬ 
fine Zll’-*- = ZllA and fi = H in the classical regime. The 
anomalous dimensions are given by: ryllA = — ^ 
and ? 7 ^ = — J^lnfl also outside the critical region, i.e., 


in general they can be scale-dependent 77 IIA = ^nd 
. Fixed points of the RG equations are identified 
by identifying constant values for all ( 77 -*-,r;!!, 77 ^). 

In the classical regime the RG equations for the 
squared field expectation value and the couplings are 
given by: 




dv 

= (2 - d - r]-^)v^ + TASdin 


ds 

dXi^2 

ds 

5^2,3 

ds 


(1 -I- v^Ai^2y 


+ N-1 


= {d-A + 277-^)Ai, 2 + Ai 27A^<i(d^) 


(1 -h 7;2 Ai,2)^ 


+ V-1 , 


— (d -I- 2z — 4 -|- 277'*‘)A2,3 -I- Ai_2A2,3T'A^(i(d'' 


9 


(1 + 7;2Ai,2)^ 


+ N-1 


The anomalous dimensions associated to the renormalization factors read 


= TA5d{0) 


' 1,2 


{i+vFx,^2y^ 


rj^=Sdiv^)^{l + 


32 




II ^ _ 5d{r]^) Ta 1_ 

2 2 v'^ kv^Xi^2 


4 -I- 7;2Ai,2(2 -I- kv'^\i^2) 

I 28 PX 12 


2 ( ’ 


1 - 


1 


1 - (k7;^Ai,2) 


(1+«2Ai,2)2 


4 -I- v^Xi^2(2 + kv'^Xi,2) 


3 ( > 


(118) 

(119) 

( 120 ) 

( 121 ) 

( 122 ) 

(123) 


and are derived for k > 0 and H > 0. Apart from a possibly nonvanishing anomalous dimension 77 ^, and the lifting 
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s = ln{k/A ) 


FIG. 7. Upper panel: Renormalization group flow for the 
parameter k evaluated at the Wilson-Fisher FP (Ta = TA,cr) 
for the = 1 scalar model. We map the infinite range of 
values K £ [0, oo) to a compact interval k i— R{1 + in 
order to illustrate the R ^ oo limit. The solution at it = 0 
(dissipative dynamics) is IR attractive as indicated by the 
arrows on the horizontal axis. The limit it —i oo corresponds 
to the relativistic fixed point and is IR repulsive. Lower panel: 
Dynamic crossover between the two scaling solutions it = 0 
and it —i oo shown for different temperatures, Ta < TA,cr 
{TA,cr- solid, blue curve). As the temperature increases the 
crossover scale is shifted towards the UV. 


FIG. 8^ The scale dependence of the renormalization param¬ 
eters and for the N = 1 scalar model is shown for 
different values of the temperature: Ta < TX.cr (TX.cr: solid, 
blue curve). As the temperature is increased towards its crit¬ 
ical value, the parameter fl displays scaling behavior, i.e., 
D ~ e”'”’ and = const., while vanishes at some finite 
scale before the asymptotic scaling regime is reached. 


To close the set of flow equations (1181 - (1251 we 
provide the flow equation for the coefficient k 


of the degeneracy of we observe two important 

differ enc es to the RG flow equations at T = 0, cf. Eqs. 
(1071 - (112): (1) Instead of H = d -I- 1 we see that 
the spatial dimension d appears in the canonical scaling 
contributions to the RG flow. (2) The loop contributions 
are not only controlled by the couplings A 1,2 and A 2,3 but 
also by the temperature. 

Note, that both anomalous dimensions 77 ^ and de¬ 
pend on the parameter 


K = 


(124) 


Its significance is clear from the discussion at the end 
of Sec. |IVB[ It defines a ratio between the parameters, 
Z^!Z} and that determine the low-energy dynam¬ 

ics of our model. In particular, we note the two inter¬ 
esting limiting scenarios: 1 ) it —f 00 : Z^^ = Z-^ = 1 , 
TaLI -t O’*' corresponding to the microscopic model with 
relativistic dynamics and vanishing fl (Ta yf 0 ), and 2 ) 
K = 0: ^11 = 0 , Z-^ > 1 , TaCi > 0 in the classical regime. 


Ok 

ds 


2(1-I- 77 ^) — 77 II — Ty*- K = (z + r]^ — (125) 


and use the relation z = 2 — 77 -*- -|- 77 ^ that defines the 
dynamic critical exponent. 

Even without solving the full set of equations, we may 
infer the dynamics of this model and determine the as¬ 
sociated fixed points. However, this discussion is some¬ 
what subtle - substituting the given expressions for the 
anomalous dimensions into Eq. (125) one might conclude 


that there is only one solution for nonnegative k, which 
is reached in the k —>■ 00 limit (see Fig. [^. Further fixed 
points exist, but all of them lie in the region of it < Oj^ 
Such negative solutions are clearly unphysical and are 


The presence of such solutions within our truncation might be re¬ 
lated to the neglect of higher order 0{lu^) terms in the frequency 
expansion of A) that may stabilize the large-frequency behav¬ 
ior of the renormalization group flow. 
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FIG. 9. = 1 scalar model {d = 3, Ta > 0): As we tune 

the temperature Ta through its critical value Ta,ci at fixed 
values of Va = Va/K = 0.078 and Aa = Aa/A = 1, we observe 
a phase transition, where the theory exhibits scaling. We 
show the scale dependence of the dynamic scaling exponent 
z = 2 — rp and find 2 ~ 2.025 in the scaling regime. 


not allowed within our truncation (which requires that 
> 0, and Cl > 0). Hence, when the renormaliza¬ 
tion group flow reaches S = 0 we stop the evolution of Z'^ 
which is then set to zero for the remaining flow. Thus, it 
appears that within our truncation only two independent 
scaling solutions exist, corresponding to k = 0 {Z'^ = 0) 
and it —)■ oo > 0). Thermal fluctuations drive the 
system towards the IR stable FP at it = 0 and we are 
able to identify a dynamic crossover (cf. Figs. mill)- This 
scenario is of course in line with our expectation that 
the dynamics is governed by nonrelativistic relaxational 
modes well below the thermal scale. We provide further 
interpretation of these results in the concluding section 
(Sec. VIIII where we also comment on how different trun¬ 
cations might affect this picture of the dynamic crossover 
between the it = 0 and it —>■ oo scaling solutions. 


Let us focus on this asymptotic scaling regime where 
the it = 0 solution takes over [Z'^ = 0). We observe that 
in this regime, Eq. (122) is identical to the anomalous 
dimension rf' of the relaxation coefficient derived in the 
context of the functional RG for Model A [5S1 [H] . We 
identify the following scaling solutions in this regime: 


(I) Gaussian fixed point: 


Ai,2 — A2,3 — 0 (126) 

and the anomalous dimensions vanish, i.e., rj-^ = 
= 0. The dynamic critical exponent is z = 2. 

(II) Wilson-Fisher fixed point (2 < d < 4): The squared 
field expectation value and the quartic coupling 
Ai ^2 show a temperature dependence: 


d\TA) = TAd\l), 

(127) 

Ai,2(Ia)=Ta^Ai,2(1), 

(128) 

A2,3(2a) = 0 , 

(129) 


at the fixed point. The anomalous dimensions are 
nonvanishing: 0. 


Although both FPs depend on the temperature, the scal¬ 
ing exponents do not. This i s easi ly ch ecked, if we sub¬ 
stitute the scaling form (127) and (128) in our result for 
the anomalous dimension 


= <5d(0) 


i^^(i)%(i) 

1 + F2(1)Ai.2(1) 


(130) 


and the dynamic critical exponent 


z = 2 - 77 -^ 


^ 2 ( 1 ) 


1 + 


1 -b 2(1) 


2 -b u2(1)Ai_2(1) 


(131) 


We therefore conclude, that as a function of temperature 
both FPs describe a continuous line of phase transitions 
which lie in the same universality class. Depending on 
spatial dimension either one of the two FPs is IR stable 
and characterizes the scaling properties at the transition. 
In particular, it is the Wilson-Fisher FP that is stable 
below the upper critical dimension d^r = 4 and we find 
that the dynamic scaling is determined by the universal¬ 
ity class of Model A (within the limits of our truncation). 

We proceed to discuss the behavior of the theory as we 
vary the temperature Ta through its critical value TA^cr 
(in units of the cutoff scale A). Our results are illustrated 


at the example of the N = 1 scalar model in d = 3 dimen¬ 
sions, cf. Figs. 0111 and As we have already discussed 
the static properties of the three-dimensional Ising uni¬ 
versality class in the previous section, we focus only the 
dynamic properties. As we have already argued, we ex¬ 
pect a dynamic crossover well below the thermal scale. 
This is verified by examining Fig. where we see that 
the renormalization coefficient Z^^ runs to zero before the 
system enters the scaling region. The scaling behavior in 
this asymptotic regime fc <C A is fully characterized by 
the Wilson-Fisher FP, for which we obtain the dynamic 
critical exponent z ~ 2.025, in the dynamic universal- 
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ity class of Model A. Our result agrees well with Monte 
Carlo estimates: z = 2.032(4) [l08], z = 2.055(10) [109], 
field-theory methods 0 = 2.0237(55) |110| . and is consis¬ 
tent with the functional RG results from Ref. [55] to the 
given order of our truncation. 

The following picture emerges: In the classical regime, 
we obtain a continuous line of phase transitions that lie 
in the d-dimensional static universality class of the 0{N) 
model. The critical exponents at the phase transition de¬ 
pend on nature of the IR stable fixed point. For N > 2 
and 2 < d < 4 the Wilson-Fisher FP is the stable one 
and we observe nontrivial scaling, i.e., t] ^ 0, z > 2. 
For d > 4 the Gaussian FP takes over and we obtain 
mean-field exponents, rj = 0, z = 2. For N = 1 (Ising 
universality class), the lower critical dimension is d = 1 
and we observe nontrivial scaling for any dimension be¬ 
tween 1 < d < 4. The dynamic critical exponent is z > 2 
for 1 < d < 4 (which is consistent with the results from 
Ref. mu and those compiled in Ref. El)- 


A. Expansion around the upper critical dimension 

As we have already noted in the previous section, the 
coupling A 2,3 vanishes exactly at the Wilson-Fisher fixed 
point. Thus, the only relevant parameters that character¬ 
ize this scaling solution are the coupling A 1.2 and the field 
expectation value v. For the following discussion it will 
prove useful to introduce the parametrization: x = v^Xi ^2 
and y = i5d(0)TAAi,2, for which fixed point equations take 
a somewhat simpler form 


0 = (2 — d — 'q'^)x -\- 


5M 


0 = (d — 4 -I- 2?7"*‘) -|- 


<5d(0) 


(1 -bx) 


2 +N -l\y, 


(132) 


9 

(1 -b x)3 


-bA^- 1 


y- 

(133) 


VII. LIMITING BEHAVIOR OF 
RENORMALIZATION GROUP EQUATIONS IN 
THE CLASSICAL REGIME 

Here, we examine the properties of the 0{N) Wilson- 
Fisher scaling solution in the classical regime (T ^ 0), 
by considering an expansion around the upper critical 
dimension of the model dcr = 4 and the behavior in 
the large-A limit. Although we do not expect to repro¬ 
duce the known results from the e- or the 1 /A-expansion 
(which is due to threshold effects in the nonperturba- 
tive functional RG), this analysis nevertheless provides 
important insights into the quality of our truncation in 
these limiting cases. 


The anomalous dimension is given by 


_L _ xy 

(l + x)2’ 

and dynamic critical exponent takes the form 


(134) 


rf(d^) 2(x-bl) _L 

5M (x + 2)2^ 


(135) 


Equations (132) and (1331 imply that the leading contri¬ 
bution to the parameters x and y in an expansion around 
the upper critical dimension e = 4 — d is of order e. Up 
to second order, we find 


X = 


y = 


N + 2 (A-b2)2(A-b33) 2 

^ ' c 


2(A-b8) 4(A-b8)3 

1 25(A + 2)__2 , ,^,,3 


0{e^ 


N + 8 2(A-b8)3 

Thus, within our truncation the anomalous dimension is given by 


Oie^). 


= 


1 A-b2 


A-b2 


2 (A-b 8)2 (A-b8)2 

to third order in e, while the dynamic critical exponent reads 


14A -b 37 

(A+ 8)2 


C>(e^ 


z = 2- 


l+Oie^) 


(136) 

(137) 

(138) 

(139) 


to the same order as the anomalous dimension. Let us comment to what extent the limiting behavior (138) and 
(1391 reproduces the known e-expansion result. For the anomalous scaling exponent ry*- observe that the O(e^) 
contribution is exactly reproduced, while higher-order terms appear to be inconsistent |115[ 1116] . This is in contrast 
to other critical exponents. In particular, from the eigenvalues of the stability matrix at the Wilson-Fisher FP, we 
determine the correlation length exponent 


v = - 
2 


1 (A + 2) (A + 2) 


4(A + 8) 8(A + 8)3 


A^ + -(152A + 574) 


+ 0(e'), 


(140) 
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Regulator function 

(Ir — -X'r^;l/^;2 

Cr — 1 

ExDonential cutoff |1()41I112| 

r-exp = (e^ - 1)“^ 

1/2 

6 In (4/3) - 1 

Litim cutoff |77| 

»'opt = il/y-l)e{l-y) 

1/2 

1/2 

Sharp cutoff |113||114| 

'^sharp ~ 1) 1 

CKi 

-1 

e-exoansion |i631111511116) 


1/2 

6 In (4/3) - 1 


TABLE 1. Numerical coefficients and Cr that determine the quality of a given cutoff in the local interaction approximation 
are compared with the 0{^) result from the e-expansion, i.e., rj^ = S'lid z = 2 + CrT]^■ Overall, the exponential 

regulator provides an 0{^) improvement for both t]^ and rf'. This is in contrast to the optimized Litim and sharp cutoff 
functions that fail to reproduce the known results. Note, that the sharp cutoff limit and its shortcomings are well-documented 
in the literature (see, e.g., [63l[77ll87llTT7]i. 


and Wegner’s exponent w 


where the coefficients and are given by 


25(iV-h2) 2 

2(nTW^ 




(141) 


which governs corrections to asymptotic scaling. Since 
there are only two independent exponents at the FP, we 
may determine the remaining exponents by use of scaling 
relations below the upper critical dimension m, e.g., 
we may apply Fisher’s scaling relation 7 = (2 — to 
determine the susceptibility exponent 


7 = 1 + 


X 


(N + 2) _ (IV+ 2) 
2(iV + 8)^'^ 4(iV + 8)3 

iV2 + i(i497V + 550) 

O 


e^ + 0{e^). (142) 


We observe that the given exponents and scaling rela¬ 
tion match those of the e-expansion m perfectly to 
linear order however, discrepancies appear for contribu¬ 
tions of order e^. In fact, linear contributions in e are 
in general independent of the chosen regulator function 
[11711120] , while this does not hold true for higher order 
contributions (at least within the employed truncation of 
the effective action). 

We demonstrate the regulator dependence of the two- 
loop result explicitly for the anomalous dimension which 
is seemingly consistent at 0{e^). Consider a momentum- 
dependent regulator R{p) = Z-^p^r[p'^/k"^) that is gen¬ 
eral arbitrary up to the specification of r = rijj), y = 
p'^/k'^. This function that implements the IR cutoff, 
should be nonnegative r{y) > 0 and we impose the fol¬ 
lowing limiting properties: limj,_>o r(j/) ~ 1 /j/ —>• 00 and 
decay sufficiently fast at infinity , i.e ., limy_).oo'r(j/) —)■ 0 
and limy_>oo ?’^(y) —t 0 (cf. Sec. |III[ ). Using a regulator 
function that satisfies these conditions but is otherwise 
arbitrary, we derive the general expression for the anoma¬ 
lous dimension to leading order 0{e^) in the local inter¬ 
action approximation of the 0{N) model: 


ri 


L 


{N + 2)e‘^ XrYr-i , 3. 

(fv + 8)2 


(143) 






I “‘'(It+W) ’ 

r ^ 2-ud 1 
Jo dy{l + r{y)) 


(144) 

(145) 


Their numeric values depend on the specific choice of the 
regulator function. Let us consider a few examples: We 
start with the optimized Litim regulator 


roptiy) = (1/y- 1)0(1 - y). 


(146) 


employed in this work, cf. Secs. |V|- |VI[ The correspond¬ 
ing coefficients read Xopt = 1 and Yopt-n = ?^/2. Thus 
we see that the ratio 2 fopt 2 ^t;i (Fopt; 2 )~^ = 1 / 2 , and we 
confirm our result Eq. ( 138[ ) for the anomalous dimen¬ 
sion. To check how (146 1 compares to other regulators, 
we consider the exponential cutoff 

?’exp(l/) = (e^ - 1)"\ (147) 


which is also frequently used and serves as a useful bench¬ 
mark. This function yields Xexp = 1/2 and the low¬ 
est coefficients l^xp;n in the series n = 1,2,3,... are 
1, 1, 31n(4/3), ..., etc. Thus, we obtain the ratio 
-^explexp;! (Uexp; 2 )~^ = 1/2 and the same leading-order 
result in the e-expansion for as for the Litim regula¬ 
tor. From these examples we conclude that although the 
result to two-loop order, Eq. (1431 is regulator dependent 
in general, that it is nevertheless possible to find a class 
of regulator functions with improved behavior at O(e^). 
A counterexample that fails to produce the known re¬ 
sult for 77 -*- to leading order is the sharp cutoff [11311114) , 
which is known to introduce strong regulator artifacts 

El [771 EH]. 


For the dynamic critical exponent z, the result from 
the e-expansion is expressed in the form 2 = 2 + cry*- 
[TT] . with the coefficient c = 0.7261 (l — 1.69e + 0(e^)) 
m- Here, we observe a clear difference to order e^, 
where we find c = 1/2 + 0 (e^) within our truncation, 
and employing the Litim cutoff. While the anomalous 
dimension ry*- jg correctly reproduced to order for this 
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choice of the regulator function, this must not necessarily 
hold for T]^. To parametrize the regulator dependence of 
our results, we write r]^ as 

= + (148) 


where the coefficients Xr and Yr-n are defined in Eqs. 
(1441 - (145). From the definition of the dynamic critical 
exponent z = 2 — r]^+ri^ = 2 + crj'^, we may therefore 
express c in the form: 


c=^-l. (149) 

Using the previously derived values, we obtain to lead¬ 
ing order in e: Copt = 1/2 for the optimized Litim cutoff, 
while for the sharp cutoff we find: Cgharp = ~1- Consid¬ 
ering this discrepancy, it is quite surprising to note that 
the same coefficient takes a different value for the expo¬ 
nential cutoff: Cexp = 61n(|) — 1 ~ 0.72609, in perfect 
agreement with the e-expansion. 

We may conclude that in general the scaling expo¬ 
nents, both static and dynamic, are regulator dependent 
at O(e^), while the result to linear order in e is correctly 
reproduced independent of the choice of cutoff (see also 
Refs. |117LI120] b Nevertheless, certain regulators might 
show improved behavior at 0{e^) (at least within the 
scope of the employed truncation of the scale-dependent 
effective action). A summary of these results is shown in 
Tab.Hl 


B. Large-A expansion 


In the same representation of the parameters (132) and 


(1331, we may consider the large-A behavior as an alter¬ 


native approach to understand the properties of our RG 
equations. In this limit, we find that 


4 - / 1 


y-{4-d)^+o(-^j, 


(150) 

(151) 


up to order 1 /N. The result for the anomalous dimension 
reads 


, (d-2)U-d)‘^ 1 

^ 1 - 


N 




while the dynamic critical exponent is given by 


z = 2- 


4(d-2) 

^2 


O 


(152) 


(153) 


In d = 3 dimensions results are readily available in the lit¬ 
erature: In the large-A limit the value for is given by 
ry*- = ^ A |12I] . Clearl y, th is result does not match the 
limiting behavior in Eq. (152) which is known to exhibit a 


regulator dependence within our finite truncation of the 
vertex expansion |105) . The difference can be explained 
by the infinite resummation of diagrams in the large-A 
result, but which cannot be accounted for in the context 
of our ansatz for the CTP effective action. Our result for 
the dynamic critical exponent should be compared with 
predictions for the dynamic critical behavior of Model A, 
for which the dynamic exponent reads z = 2 + crj -^, with 
c = 0.5 in the large-A limit [TTl I 122 j . 


VIII. QUANTUM-CLASSICAL TRANSITION 


Although in principle the flow equations can be 
derived also for the intermediate temperature region 
(0 < Ta ;< 1 ) we find that they take a rather compli¬ 
cated form and therefore require a more elaborate nu¬ 
merical treatment. We postpone this to future work. 
Here, we examine possible scenarios of the quantum-to- 
classical crossover based on our understanding of the low 
and high-temperature limits, developed in Secs. and 
|VI[ The general picture of the phase structure is the fol¬ 
lowing: The d-dimensional 0{N) theory is characterized 
by a continuous line of classical phase transitions that 
terminates at T = 0. Along this line the dynamic scal¬ 
ing exponent \s z > 2. At the end point, the system 
lies in a different universality class, namely that of the 
(d -|- I)-dimensional Euclidean theory with z = 1. 

Essential to the picture of the dynamic crossover is the 
observation that both in the quantum and classical case 
we may find two fixed points: the Gaussian FP (ry = 0) 
and Wilson-Fisher FP (?y 7 ^ 0). The existence of a non¬ 
trivial IR stable FP depends subtly on the dimensionality 
and temperature of the system. In contrast, the dynamic 
scaling properties turn out to be much more simple - 
they completely decouple from the statics for any T > 0. 
This is immediately clear from the inspection of the RG 
flow Eqs. (118) - (125) where we see that the parameter 
R and the anomalous dimensions 77 ^ and 77 H do not feed 
back into the RG flow of the static quantities. In fact, we 
have used this property explicitly to determine the dy¬ 
namic FPs of the RG flow, i.e., R = 0 and k = 00 for any 
T > 0. Thus, the problem of determining the dynamics 
at nonvanishing temperature reduces to the problem of 
understanding how the Wilson-Fisher FP and in partic¬ 
ular how rj^ depends on the effective dimension d^ (cf. 
Sec. IVG 2). Between the quantum and classical regime, 
Ty-*- varies continuously with dg, from d-|-1 to d. This fact 
can be understood from the perspective of Euclidean fi¬ 
nite size systems where one observes a continuous dimen¬ 
sional crossover |123j . 

Possible scenarios for the RG flow are illustrated in 
Fig. at the example of the A = 1 scalar theory. 
We examine the scale-dependent anomalous dimension 
Ty*- = —(d/ds)lnZ-^ which provides information on rel¬ 
evant fluctuations in different regions of the RG flow. 
The maximum achievable correlation length ^ 
is determined by how close the microscopic parameters 
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FIG. 10. = 1 scalar model for spatial dimensions d — 2 and d = 3: Schematic scale dependence of the anomalous dimension 

. We tune the bare field expectation value Ua as a function of the temperature, so that the theory exhibits scaling in the 
fc —> 0 limit. We assume that the bare coupling Aa is hxed at a given value. In d = 2 dimensions, where the theory admits a 
QCP characterized by nontrivial scaling the system features a crossover scale to the IR scaling region (left panel). This scale 
is set by the quartic coupling. As we turn on temperature, quantum critical fluctuations die out and thermal fluctuations 
take over. In d = 3 dimensions (right panel) the QCP lies in the (3 + l)-dimensional (Euclidean) universality class which is 
characterized by mean-field scaling - fluctuations become important only at nonvanishing temperatures. The curves are shown 
to scale with = 1/4 (d = 2) and = 0.03639(15) (d = 3) for the Ising model. 


have been tuned to their respective critical values, and 
we assume that ^ oo in the limit fc —?> 0. For small 
temperatures (0 < Ta ^ 1), the RG flow enters a regime 
where quantum fluctuations are present but do not play 
out their full effect. Instead, close to the classical phase 
transition the characteristic scale at which fluctuations 
become important are fully determined by the temper¬ 
ature. We point out that the possible situations illus¬ 
trated in Fig.[^have to be modified in the case of a sys¬ 
tem with continuous 0{N) symmetry, where for d = 2 
and N = 2 the system experiences Kosterlitz-Thouless 
transition, while a thermal phase transition is absent for 
N>3. 


CONCLUSIONS 

In this work, we have applied the nonperturbative 
functional RG in the real-time formalism to the prob¬ 
lem of the low-energy dynamics of a relativistic 0{N) 
theory. Our approach relies on a truncation of the scale- 
dependent GTP effective action and we have shown that 
our ansatz admits two possible fixed points at T > 0 
that determine the low-energy dynamics of the order pa¬ 
rameter. Only one of them is IR stable which charac¬ 
terizes the relaxational dynamics of a dissipative mode. 
We find that the dynamic universality class associated to 
this scaling solution is identical to that of Model A, for 
which z = 2 + cr]-^ [II1IIQ3]. This observation is based 
on an exact equivalence of the flow equations derived in 
this work in the IR scaling regime, with those obtained 
in the framework of the effective model given in Refs. [28] 
and |29|. At T = 0, we find only one fixed point, with 
unitary dynamics z = 1 protected by Lorentz symme¬ 
try. We show that there is a continuous crossover that 


smoothly connects the microscopic theory {z = I) with a 
nonrelativistic effective theory (z > 2) that corresponds 
to Model A, if the theory is tuned to criticality. 

We emphasize that the existence of the Model A fixed 
point in the high-temperature limit (within our trun¬ 
cation of the scale-dependent effective action) is essen¬ 
tially required by consistency of our ansatz (cf. Sec. jVIj ). 
The stability of the theory (i.e., > 0, Q > 0) 

had to be imposed to rule out a renormalization group 
flow to negative values of that are clearly unphysi¬ 
cal. Although the relativistic UV fixed point seems to 
be unstable, which is certainly expected at nonvanish¬ 
ing temperatures, the apparent lack of a solution where 
111 = 0 is reached asymptotically, is certainly surpris¬ 
ing. We suspect that this is an artifact of our truncation 
that includes only a linear contribution in the frequency- 
expansion of the self-energy and is therefore insufficient 
to capture the large-frequency behavior of correlation or 
response functions. To reach a conclusion on this issue 
other types of approximations will need to be employed 
that take into account the full frequency-dependence of 
the self-energy and vertex functions |124H127] . 

Within our truncation, we were able to observe a dy¬ 
namic crossover induced by thermal fluctuations. How¬ 
ever, as far as the low-energy limit is concerned, one 
might also also imagine a different scenario, where the 
system features a crossover to the nonrelativistic limit 
well above the thermal scale. This would have dramatic 
consequences, since the corresponding nonrelativistic ac¬ 
tion conserves particle number |I19L 1128] . Thus, the 
question that remains is: What is the nature of the non¬ 
relativistic limit? We see that an answer to this ques¬ 
tion is intimately tied to that of the dynamic universal¬ 
ity class for the given microscopic theory. In either one 
of the above scenarios the nonrelativistic action will fea- 
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ture different symmetries, that determine the presence 
or absence of further slow modes at the phase transition. 
To capture these possibilities different truncations of the 
scale-dependent effective action will have to be consid¬ 
ered, that are able to take into account conservation laws 
and Ward identities consistently. 

Thus, within our current ansatz for the CTP effec¬ 
tive action Eq. (60) it therefore not possible to pro¬ 
vide a definitive answer to the question of the dynamic 
universality class of the 0{N) theory. It is only pos¬ 
sible to follow the transition from the microscopic co¬ 
herent dynamics (without dissipation) to the situation 
where the dynamics is dissipative nonrelativistic, with¬ 
out additional conserved quantities - which corresponds 
to Model A dynamics in the critical region. The simplest 
possible extension of this work would be to allow for a 
diffusive relaxation of the order parameter in addition 
to a purely dissipative behavior. This might allow one 
to distinguish between two dynamic fixed points associ¬ 
ated to the universality classes of Model A and Model B 
in the Halperin-Hohenberg classification. A correspond¬ 
ing ansatz for the CTP effective action should include 
a momentum-dependent as well as a field-dependent ki¬ 
netic coefficient, i.e., flab = ^ab{(l>) + ^'abWP^ + 

We do not expect this to be sufficient however, to address 
the dynamic properties of the 0{N) model in the low- 


energy limit. Certainly, the interplay of massless Nambu- 
Goldstone modes, and the possible presence of particle- 
number and energy-conservation allows for very compli¬ 
cated dynamics with competing mode couplings. To dis¬ 
tinguish between generic dynamic universality classes, we 
expect that a nonlocal expansion in the vertices (both in 
space and in time) is necessary (see, e.g.. Refs. |124P 
iMunni)- Such a treatment is in particular required to 
address the issue of the nature of the nonrelativistic limit. 
We refer to future work which is currently in progress. 
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Appendix A: Renormalization group equations 


The flow equations for the (m -|- n)-point functions T^”*’"^ are derived from Eq. (34) by taking functional deriva¬ 
tives with respect to the fields and evaluating these expressions in the homogeneous background held conhguration 
(j)a = vSai, V > 0, and (f> = 0. Within our truncation, we assume that all (m-l-n)-point functions pl’"'"), with m-|-n > 4 
are zero. Furthermore, 4-point functions are assumed to be local in both space and time, while this not true in gen¬ 
eral for 2- or 3-point functions. The presence of a n onlinear huctuation theorem implies that only two of the 4-point 
vertices and pCJ) are nonzero (cf. Sec. IVD). With these approximations we obtain the how equations for the 
(amputated) IPI two-point function pC^) 

^pd.i) (p) = _ ^ ^ Tr {F(g)r(3.i) + 2F(9)r("’i) (-P, p + q)Re G^(p + q) (p, q )} , (Al) 


and the (nonvanishing) IPI 3- and 4-point functions 


_^r(3.i) ^ _ 
ds 

_^r(i.3) ^ _ 

ds 


ds dR 


dR d 
ds dR 

ds dR 


{F(g)r(2.i)(_5^p + p' + q) ReG^(g)r(3’i) 
Tr {F(g)r(3d) Re G^(g)r(3’i) ^ | ^ 

Tr {F(g)r(34) + ... } . 




(A2) 


(A3) 


(A4) 


We dehne the four-momenta p = (p^,p), q = {q^,q)-, etc. and the frequency-momentum integration = 

( 2 T 1 -U +1 / dq^d’^Qj while the trace Tr {• • • } runs over internal held indices. In Eqs. (Al) - (A4) we have pulled out the 
derivative operator, {dR/ds){d/dR) which acts on the regulator function. This greatly simplihes the evaluation of 
the terms on the RHS. In particular, using a frequency-independent regulator (at nonvanishing temperature), we may 
calculate the frequency integral on the RHS before evaluating the i?-derivative. Finally, those terms that are omitted 
in Eqs. (A2) - (|A4[) (denoted by the ellipsis) correspond to all possible contributions (allowed by the diagrammatic 
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rules) where the 4-point vertices of the first term on the RHS are substituted by the dumbbell diagram as illustrated 
in Fig. 11 This is done explicitly in Eq. (Al). 



FIG. 11. Substitution rule illustrated at the example of the amputated 4-point function in terms of two 3-point vertices 

(cf. Eq. (All). External lines indicate the inflowing momenta and the inserted fields in the RA-representation (full 
lines correspond to (j) field insertions, while dashed lines to insertions of the 0 field). The notation refers only 

the incoming momenta. The full-dashed lines which connect the two vertices on the RHS corresponds to (nonperturbative) 
retarded/advanced propagators, i.e., 


Equations (|A2 ) 
p —■ 0 and p' 


(A4) provide the RG flow of the generalized potential in the limit of zero external momenta, i.e., 
0. Similarly, we obtain the flow equations for the wavefunction renormalization Z-^ and the dynamic 


parameters 12 and of our model by a suitable projection onto the flow equation for 


^ - 5 'ia f ^ T.- {r(,)r<«>(-p, p + ,) Re cqp + 5)} . (A6) 

^ = -i/?h)m ^^^^Tr|F(q)r(2>i)(-p,p-hg)ReG^(p + q)r(2.i)(p^g)| . (A7) 


In practice, one usually considers only the contribution from the Goldstone modes on the RHS, since it is understood 
that they give the dominant contribution to the wavefunction renormalization Z~^ [6311105j (and similarly to the 
dynamic coefficients Z'^ and H). Here, we follow the same strategy and therefore, instead of evaluating the trace we 
carry out a projection onto the Goldstone modes in Eqs. (A5) - (A7). The diagrams that we evaluate to obtain the 
RG flow of the parameters are shown in Fig. 


12 


p + q,r P + q,^ 



FIG. 12. Diagrams that contribute on the RHS of Eqs. (A5| - (A7l. Full-dashed lines corresponds to (nonperturbative) 
retarded/advanced propagators G^^^, while the full lines denote statistical propagators F. Full/dashed external legs denote 
<()/0-field insertions. Next to the momenta we indicate the relevant modes indicated (r: radial mode, g: Goldstone mode). 
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